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ABSTRACT 



We use data obtained with the Far Ultraviolet Spectroscopic Explorer (FUSE) 



to determine the interstellar abundances of DI, NI, 01, Fell, and H2 along the 
sightlines to WD 1034+001, BD+39 3226, and TD1 32709. Our main focus is on 
determining the D/H, N/H, O/H, and D/O ratios along these sightlines, with log 

1^ ' ^(H) > 20.0, that probe gas well outside of the Local Bubble. Hubble Space Tele- 

■ scope (HST) and International Ultraviolet Explorer (IUE) archival data are used 

to determine the HI column densities along the WD 1034+001 and TD1 32709 
>• ■ sightlines, respectively. For BD+39 3226, a previously published iV(HI) is used. 

We find (D/H)xl0 5 = 2.14 ± 1.17 ± and 1.86 ± and (D/O)xl0 2 
= 6.31 ± \ f 8 , 5.62 ± and 7.59 ± 2 ;^, for the WD 1034+001, BD+39 3226, 
and TD1 32709 sightlines, respectively (all la). The scatter in these three D/H 
ratios exemplifies the scatter that has been found by other authors for sightlines 
with column densities in the range 19.2 < log iV(H) < 20.7. The D/H ratio 
toward WD 1034+001 and all the D/O ratios derived here are inconsistent with 
the Local Bubble value and are some of the highest in the literature. We discuss 
the implications of our measurements for the determination of the present-epoch 
abundance of deuterium, and for the different scenarios that try to explain the 
D/H variations. We present a study of D/H as a function of the average sightline 

^ ■ gas density, using the ratios derived in this work as well as ratios from the lit- 

erature, which suggests that D/H decreases with increasing gas volume density. 
Similar behaviors by other elements such Fe and Si have been interpreted as the 
result of depletion into dust grains. 

Subject headings: ISM: Abundances — ISM: Evolution — Ultraviolet: ISM — 
Stars: Individual (WD 1034+001, BD+39 3226, TD1 32709) 
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1. INTRODUCTION 

The present day abundance ratio of deuterium to hydrogen places important constraints 
on Big Bang nucleosynthesis (BBN) and the chemical evolution of galaxies. Because deu- 
terium is only produced in appreciable amounts in primordial BBN and destroyed in stellar 
interiors, the measurement of D I/H I in the interstellar medium (ISM) places a lower limit 
on the primordial abundance of deuterium. In addition, by comparing the ISM abundance 
of deuterium to its abundance in high-redshift intergalactic gas we should be able to better 
understand the effects of astration and chemical evolution of galaxies. Measurements of the 
D/H ratio in intervening clouds of gas seen toward distant quasars have yielded a range of 
values D/H = (1.65 - 4.0) x KT 5 (O'Meara et al. 2001; Pettini & Bowen 2001; Levshakov et al. 
2002, and references therein). Kirkman et al. (2003) measured D/H = (2.42 ± alfDxlO 5 
toward Q1234+3047 (z = 2.526). These authors calculate the value they believe is the best 
estimate of the primordial D/H abundance, D/H prim = (2.78 ± o^xlO" 5 (lcr errors in the 
mean), by taking the weighted mean of five D/H measurements toward QSOs ranging from 
z = 2.079 - 3.572. This value is in good agreement with that determined from the cos- 
mic microwave background measurements performed with WMAP and previous missions 
(Spergel et al. 2003). Sembach et al. (2004) found D/H = (2.2 ± 0.7) xl0~ 5 for Complex C, 
a high velocity cloud falling into our galaxy, which has low metallicity and has presumably 
experienced more stellar processing than the gas seen toward QSOs. 

Measurements of D/H in the local ISM have been made with Copernicus (e.g. Rogerson 
& York 1973), HST, (e.g. Linsky et al. 1995), IMAPS, (Jenkins et al. 1999; Sonneborn et al. 
2000), and more recently FUSE (e. g. Moos et al. 2002, and references therein). A nearly 
constant ratio of D/H = (1.5 ± 0.1) x 10~ 5 (la on the mean) has been obtained in the Local 
Interstellar Cloud (LIC) by Linsky (1998); recent measurements inside the Local Bubble 
(LB, log iV(HI) < 19.2 Sfeir et al. 1999) appear to be consistent with a single value for D/H 
in the LB (Moos et al. 2002; Oliveira et al. 2003). From a compilation of measurements from 
the literature Wood et al. (2004) derive (D/H) LB = (1.56 ± 0.04) xl0~ 5 . Using (D/0) LB 
= (3.84 ± 0.16)xl0- 2 from Hebrard & Moos (2003) and (0/H) LB = (3.45 ± 0.19) xl0~ 4 
from Oliveira et al. (2005) one can derive indirectly (D/H) LB = (1.33 ± 0.09) xlO" 5 . The 
direct and indirect determinations of (D/H)l B are only consistent when one considers the 
2a uncertainties in both quantities. A detailed discussion of the possible causes for this 
discrepancy can be found in Hebrard & Moos (2003). 

Outside the Local Bubble, however, there still is not a consistent picture of the D/H 
behavior. Measurements performed along sightlines probing gas outside the LB suggest 
variations of the interstellar D/H ratio beyond the Local Bubble, at the distance of a few 
hundred parsecs. D/H varies for 19.2 < log iV(H I) < 20.7 (e.g. Jenkins et al. 1999; Sonneborn 
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et al. 2000; Friedman et al. 2002) and apparently remains constant for log iV(HI) > 20.7, 
albeit with a lower value. The possibility of a trend towards low D abundance at high H 
column densities was noted in the D/O survey by Hebrard & Moos (2003). Wood et al. 
(2004), using additional measurements and values in the literature, confirmed the trend 
for D/H. So far, only five measurements for sightlines with log iV(HI) > 20.7 have been 
published. Hoopes et al. (2003), Wood et al. (2004) and Hebrard et al. (2005) measured 
D/H along five extended lines of sight (d > 500 pc). The weighted mean of these five 
measurements yields D/H = (0.87 ± 0.08) xlO -5 , in clear disagreement with the LB value. 

The value of the present-epoch Milky Way abundance of deuterium, (D/H)pe, is still 
an issue subject to debate. Two opposite explanations have emerged. (Hebrard & Moos 
2003) argue for a value of the current-epoch D/H ratio lower than D/H L b, on the basis of the 
generally low values found for more distant sightlines. On the other hand, Linsky et al. (2005, 
in prep) argue that an important fraction of deuterium along many sightlines is trapped 
in a population of grain material and hence the current-epoch deuterium abundance (gas 
plus grains) is (D/H) PE > (2.19 ± 0.27) xl0~ 5 . The astration factors, (D/H) prim /(D/H)p E , 
derived from both of these scenarios (~4 and ~1.2, respectively) challenge current models 
of galactic chemical evolution. 

Because the O/H variations in the diffuse ISM from the LB out to 1 Kpc, are small, 
consistent with the measurement uncertainties (Meyer et al. 1998; Andre et al. 2003; Car- 
tledge et al. 2004; Oliveira et al. 2005), OI has been used as a proxy for HI (for a more 
detailed discussion see Hebrard & Moos 2003, and references therein). Thus D/O and D/H 
are expected to trace each other in the diffuse ISM. However, how well this works at high 
column densities is a subject of discussion (Hebrard 2006, and §7.1 of this paper). 

At the present time, definitive conclusions are prevented by the limited number of D/H 
and D/O measurements along high column density sightlines. In this work we present D/H 
and D/O measurements along three sightlines (WD 1034+001, BD+39 3226, and TD1 32709) 
with log iV(HI) > 20.0. Using data obtained with FUSE, IUE, and the Goddard High 
Resolution Spectrograph (GHRS, onboard HST) we derive column densities of atomic and 
molecular species which are then used to determine important ratios (D/H, N/H, O/H, etc.) 
that can be compared to values in the literature. 

This paper is organized as follows. The three targets used for the analyses are described 
in §2. The observations and data processing are presented in §3. In §4 we determine the 
column densities of HI, DI, NI, OI, Fell, and H 2 along the WD 1034+001, BD+39 3226, 
and TD 132709 sightlines. Lower limits to the column densities of other elements are also 
reported. The D/H, D/N, D/O, N/H, O/H, and O/N ratios for the three sightlines are pre- 
sented in §5 and discussed in §6, where evidence for above average D abundance is presented 
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for the WD 1034+001 and TD1 32709 sightlines. The D/O ratios for the three sightlines are 
substantially higher than the mean LB value, and thus they considerably increase the scatter 
in the published D/O ratios. Implications of the derived ratios for hypothetical mechanisms, 
such as dust depletion and infall of metal poor gas, that try to explain the D/H variations in 
the Galactic disk are considered in §6.4. Our derived D/H, D/O, and O/H ratios are com- 
pared with previously published ratios in §7.1. In §7.2 we consider the relationship between 
D/O and Fe/O. In §7.3 we perform a study of D/H as a function of the average sightline gas 
density (N(B.)/d). Our study indicates that D/H is constant up to densities of 0.10 cm -3 
and decreases with increasing density after this point, similarly to what has been observed 
for other elements such as Fe and Si. We summarize our findings in §8. All uncertainties are 
quoted at the la level unless stated otherwise. 

2. THE TARGETS 

The properties of the three stars are listed in Table 1. Below we discuss each star in 
detail. 



2.1. WD 1034+001 

WD 1034+001 is a hot DO white dwarf with T cff =100,000 ± ^ooo K > lo g 9 = 7 - 5 
± 0.3, at a distance of 155 ± pc (photometric, Werner et al. 1995) in the direction / 
= 247.55° and b = +47.75°. This star was first observed in the ultraviolet with IUE by 
Sion et al. (1985), who identified several photospheric features (CIV, NV, and OV) as 
well as interstellar features (NI, CII, Sill, Mgll, Sill, and Si III). These authors derived 
fPH — fisM ~50 km s -1 (vph and t>isM correspond to stellar and interstellar absorptions, 
respectively). Non-stellar absorption by CIV was also detected along this sightline by Sion 
et al. (1985), however the resolution of the data did not allow them to decide whether this 
absorption had an interstellar or circumstellar origin. 

Werner et al. (1995) used GHRS onboard HST to observe this target at a resolution 
of ~20,000. In addition to stellar absorption by Hell, CIV, NV, OV, Si IV, FeVI, and 
FeVII, they also identified absorption by highly ionized species of CIV, NV, and Si IV. 
With resolution similar to that of Sion et al. (1985), Werner et al. (1995) could not determine 
whether the absorption had an interstellar or circumstellar origin. Using H I Lya observations 
obtained with GHRS in conjunction with NLTE model atmospheres, Werner et al. (1995) 
determined log iV(HI) = 20.05 ± 0.20. In this work we reanalyze these HI Lya observations 
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together with FUSE data to place tighter constraints on iV(HI) (see §4.3.1). 

Using Sloan Digital Sky Survey data, Hewett et al. (2003) discovered a region of ionized 
gas around WD 1034+001, seen in H Q , Hp, [N II], [O III], and [S II] emission, with a diameter 
greater than 2° and identified it with a planetary nebula, Hewett 1. This is one of the largest 
planetary nebula in the sky, and also the first to be unambiguously associated with a DO 
white dwarf. Using data obtained with the Southern H- Alpha Sky Survey Atlas, Rauch et al. 
(2004) found two extended emission structures surrounding WD 1034+001, in addition to 
the planetary nebula discovered by Hewett et al. (2003). The inner halo has a linear diameter 
of 16.2 ± 45 x 24.3 ± gg pc while the wider shell has a linear diameter of 27x43 pc (Rauch 
et al. 2004), at the quoted distance for WD 1034+001 (see Table 1). During its evolution, 
WD 1034+001 is likely to have passed through a much hotter phase with T eS = 150,000K, 
roughly 30,000 years ago (Rauch et al. 2004). These authors estimate that WD 1034+001 is 
most likely the exciting star of the nebula and the halo given that the recombination time 
scale in such a low-density (n e = 3 cm~ 3 ) nebula is longer than 30,000 years. 



2.2. BD+39 3226 

BD+39 3226, first identified as a sdO by Berger & Fringant (1978), lies at a distance 
of 290 ± \f pc (Hipparcos) in the direction I = 65.00°, b = +28.77°. Bluhm et al. (1999) 
determined the abundances of the atomic (HI, DI, etc.) and molecular species (H 2 ) along 
the line of sight using data obtained with IUE and ORFEUS II, and determined vph - fisM ~ 
—255 km s -1 . However, the spectral resolution of only ~30 km s _1 for both the IUE and 
ORFEUS 7/datasets led to fairly large uncertainties on the column densities determined with 
those data. The new FUSE data with a resolution of ~20 km s" 1 and better S/N ratio allows 
us to determine more accurate column densities and in addition, offers us the possibility of 
comparing column densities determined from data obtained with different instruments. 

In this work we use the FUSE data to revisit the column densities along the BD+39 3226 
sightline for the species that have transitions in the FUSE bandpass. Since the Lya transition 
of HI is not covered by FUSE we use the value reported by Bluhm et al. (1999). These 
authors used Lyct observations from IUE and ORFEUS II, in conjunction with a stellar 
model with T eff = 45,000 K and logg = 5.5, to determine log iV(HI) = 20.08 ± 0.09. 
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2.3. TD1 32709 

TD1 32709 (UVO 0904-02) was first identified as a sdO by Berger & Fringant (1980) 
using data obtained with the Ultraviolet Sky Survey Telescope S2/68. Using optical data in 
conjunction with NLTE stellar atmosphere models, Dreizler (1993) determined T e g- = 46,500 
± 1,000 K and log g = 5.55 ± 0.15. The corresponding photometric distance is d = 520 ± 
90 pc (z = 245 ± 40 pc). IUE observations covering the Lya region in combination with 
FUSE data and stellar models allow us to determine iV(HI) along this sightline (see §4.3.2). 
From our work we determine Vph — fisM = — 13 km s" 1 (see §4.2.3). 

3. OBSERVATIONS AND DATA PROCESSING 

3.1. FUSE Observations 

The FUSE observatory consists of four coaligned prime-focus telescopes and Rowland- 
circle spectrographs that produce spectra over the wavelength range 905 - 1187 A, with a 
spectral resolution of ~15 - 20 km s _1 (wavelength dependent) for point sources. Details 
about the FUSE mission, its planning, and on-orbit performance can be found in Moos 
et al. (2000) and Sahnow et al. (2000). 

Table 2 summarizes the FUSE observations of the three targets. The data were ob- 
tained through the medium size aperture (MDRS, 4" x 20"), in histogram (H) or ttag modes 
(T). The two-dimensional FUSE spectra are reduced using the CalFUSE pipeline v2.4.1 
or V2.4.2 1 . The processing includes data screening for low quality or unreliable data, ther- 
mal drift correction, geometric distortion correction, heliocentric velocity correction, dead 
time correction, wavelength calibration, detection and removal of event bursts, background 
subtraction, and astigmatism correction. The spectra are aligned by cross-correlating the 
individual exposures over a short wavelength range that contains prominent spectral features 
and then coadded by weighting each exposure by its exposure time, using the CORRCAL 
software developed by S. Friedman. All the spectra are binned to three pixel samples, or ~20 
mA, for analysis (the line spread function, LSF, is about 11 pixels or ~70 mA wide). For 
each target, all the observations were coadded in order to increase the S/N of the dataset. 
The S/N per unbinned pixel in the SiC IB channel around 940 A is 8, 26, and 14, for 
WD 1034+001, BD+39 3226, and TD1 32709, respectively. 

Figures 1, 2, and 3, present the FUSE spectra of the three targets. The most prominent 



1 Thc CalFUSE pipeline reference guide is available at http://fuse.pha.jhu.edu/analysis/pipelinc_refcrcnce.litml 
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interstellar and photospheric lines are labeled; photospheric lines are indicated by [], and 
H 2 absorption by dashed vertical lines. 

4. ANALYSIS 

Whenever possible we use apparent optical depth, curve of growth and profile fitting 
methods (hereafter AOD, COG, and PF, respectively) to determine column densities (see 
e.g. Oliveira et al. 2003, and references therein for a further discussion of these methods). 
AOD is used on weak transitions to determine column densities of species that have multiple 
transitions in the FUSE bandpass, which are then compared to the values derived with the 
other methods. For species that only have saturated transitions in the FUSE bandpass we 
use the AOD technique to place lower limits on the column densities of those species. 

We use the COG method to determine the column densities of species that have several 
transitions in the FUSE bandpass (one or more of these transitions must be unsaturated). We 
measure the equivalent widths of all the non-blended transitions, in all the FUSE channels 
where those transitions are covered, and compare them to check for inconsistencies (due to 
fixed pattern noise, for instance). A single component Gaussian curve-of-growth is then fit 
to the measured equivalent widths of the atomic species, allowing us to determine N and b, 
for each species. For H 2 we fit a single-component Gaussian COG to the equivalent widths 
of all the J levels, simultaneously. We determine then N(J) and b, common to all J levels. 

With PF we fit a single absorption component to one or more non-saturated transitions 
of the species for which N is being sought. We use the profile fitting code Owens, developed 
by Martin Lemoine and the French FUSE Team. Details of how the PF method is used for 
the analysis of each sightline studied in this paper are given below. More information about 
the fitting code can be found in Lemoine et al. (2002) and Hebrard et al. (2002). 

We follow the procedures outlined in Oliveira et al. (2003) to use these techniques and 
determine the uncertainties in N associated with each method. We introduce a modification 
from Oliveira et al. (2003) regarding the PF technique. A single-Gaussian with a FWHM of 
10.5 pixels (constant across all wavelengths and channels and during the fitting procedure) is 
used here to describe the FUSE line spread function (LSF), whereas previously the LSF was 
a free parameter of the fitting procedure during the initial stages of the fit (see Oliveira et al. 
2003, for more details). However, since several studies have shown that column densities 
obtained with fits using variable and fixed LSFs are similar (e.g. Williger et al. 2005) we 
prefer to decrease the number of degrees of freedom of the profile fitting routine by using a 
fixed LSF. 
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All the atomic lines used in the AOD, COG, and PF analyses are shown in Table 3 along 
with log fX for each transition. We use the compilation of Morton (2003) for the atomic data 
and the ones by Abgrall et al. (1993a) and Abgrall et al. (1993b) for the H 2 data. For each 
star, A, C, and P denote transitions that are used with the apparent optical depth, curve 
of growth, and profile fitting methods, respectively. Table 4 presents the equivalent widths 
of the atomic lines used with the COG method for each sightline. After column densities 
and uncertainties are determined with these different techniques the results are examined to 
check the consistency between the different methods. The adopted results are a subjective 
compromise between the values obtained with the several methods used; the uncertainties 
adopted are such that they in general include the extreme values obtained with the different 
methods. No attempt was made to combine the different results in any statistical way. The 
results of the individual methods are also presented (see § 4.2). 

In § 4.1 we determine the H 2 column densities along the three lines of sight, atomic 
column densities are discussed in § 4.2. Determination of the HI column density along the 
WD 1034+001 and TD1 32709 sightlines is discussed in §4.3. 

4.1. Molecular Hydrogen 

4.1.1. WD 1034+001 

Along this sightline we detect absorption arising from the H 2 rotational levels J = 
- 3. We use AOD, single-6 COG and PF to derive the column densities of the J = - 3 
levels. From the single-6 COG fit we derive &cog — 3.6 ± q| km s _1 . Profile fitting of the 
different J levels is done simultaneously with the atomic species along this sightline, under 
the assumption that the molecular gas does not trace precisely the atomic gas (i.e. atomic 
and molecular species are fit in different absorption components). Although the velocities of 
the atomic and molecular component are a free parameter of the fit, no significant velocity 
differences are found between the two components. The adopted column densities for each 
J-level are summarized in Table 5. We adopt log iV(H 2 ) = 15.72 ± [j jj. The excitation 
diagram for H 2 along this sightline is presented in the top panel of Figure 4. We derive T 02 
= 341 ± 75 K. 



4.1.2. BD+393226 

For this sightline we detect absorption arising from the H 2 rotational levels J = - 5. 
We use AOD, single-6 COG and PF to derive the column densities of the J = - 3 levels. 
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The J = 4, 5 levels have very low column densities and we use only PF to determine the 
column densities associated with these two J levels. 

From the single-6 COG fit to the J = - 3 levels we derive b C oc — 4.3 ± ^ km s _1 . 
Similarly to WD 1034+001, profile fitting of the J = - 5 levels is done simultaneously with 
the atomic species, but in different absorption components. Table 5 presents the adopted 
column densities for this sightline. We adopt log iV(H 2 ) = 15.65 ± qo7- The excitation 
diagram for H 2 along this sightline is presented in the middle panel of Figure 4. The popula- 
tions of the levels with J > 2 are larger than expected if the different levels were populated 
according to a Boltzmann distribution with an excitation temperature corresponding to Tqi 
(indicated by a dashed line in Figure 4). This nonthermal excitation is the result of pump- 
ing by UV photons followed by cascading transitions down through the various rotational 
levels (Black & Dalgarno 1976). The high J lines can typically be fit by a single excitation 
temperature. However, in the case of this sightline a single temperature fit to J = 2 - 5 
clearly overestimates the column densities of the J = 3 and J = 4 levels. We find that the 
distribution of H2 through the different J levels is better described by several temperatures. 
We derive then T 01 = 104 ± 27 K, T 13 = 200 ± 15 K, and T 35 = 953 ± 158 K. 

4.1.3. TD 132709 

We use AOD and PF to determine the column densities of the H 2 rotational levels J 
= 0-3. For this sightline, the low column density of each J level implies that the number 
of absorption lines that could be used with the COG method is small, and hence we do not 
perform a COG fit for H 2 along this sightline. The consistency of PF results is checked against 
that of the AOD method. Table 5 presents the adopted N(J) for this sightline. We derive 
iV(H 2 ) = 14.48 ± oil- The bottom panel of Figure 4 presents the H 2 excitation diagram for 
this sightline. All the J levels are consistent with a single temperature corresponding to T 02 
= 292 ± 83 K. 

No HD is detected along any of these sightlines. Typical HD /H 2 ratios along these types 
of diffuse sightlines are of the order of lxl0~ 5 or less. Considering iV(H 2 ) quoted in Table 
5 we expect log JV(HD) < 11. 

4.2. Atomic Species 

In this section we discuss the determination of the atomic column densities along the 
three sightlines. However, determination of iV(HI) along the WD 1034+001 and TD1 32709 
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sightlines will be discussed separately in §4.3. Tables 6, 7, and 8 present the column densities 
for D I, N I, O I, and Fe II obtained with the different methods discussed above. The adopted 
column densities for the atomic species along the three sightlines are summarized in Table 
9. 



4.2.1. WD 1034+ 001 

We use the AOD method to derive lower limits to the column densities of CII*, C III, 
Nil, Sill, PII, and Arl, because all the transitions of these species in the FUSE bandpass 
are saturated. The AOD method is also used for weak transitions of D I, N I, O I, and Fe II, 
to determine column densities which are then compared to the column densities determined 
with the other methods. 

We fit a single-Gaussian COG to the measured equivalent widths of D I, N I, O I, and 
Fell. The four COGs yield 6-values in the range 6 = 5-6 km s _1 . Figure 5 presents the 
curve of growth for O I along this sightline. 

With PF we fit a single absorption component to non-saturated lines of DI, NI, 01, 
and Fe II. In order to model the continuum in the vicinity of the D I lines an extra absorption 
component with HI only, is also included in the fit. This component is not used to derive 
iV(HI) along this sightline. No stellar model is used in the fitting process because the 
photospheric spectrum is generally flat in the regions of interest. As discussed above, H 2 is 
also included in the fit, in a different absorption component. The DI lines used in this 
analysis have a minor blending with H 2 , hence iV(D I) is not affected by the assumed iV(H 2 ). 

Along this sightline we also detect absorption by OVI (A1032, 1038), S III (A1012), 
and SIV (A1063). OVI absorption has two components, separated by ~ 50 km s _1 . The 
strongest component falls at the stellar velocity derived by Sion et al. (1985) and is likely 
of stellar origin. Absorption by S III, S IV, and the weaker O VI absorption component fall 
at the velocity of the atomic gas (DI, NI, 01, and Fell) along this sightline. However, we 
cannot conclude with certainty whether these features have an interstellar or circumstellar 
origin. Note that as discussed in §2.1, two large ionized halos have been discovered in the 
vicinity of WD 1034+001. The O VI lines are weak and blended with stellar O VI absorption; 
we derive N(0 VI) with PF, following the method described in §4. The S III absorption line 
is strong and likely saturated. We derive only a lower limit to its column density using the 
AOD technique. The S IV line is blended with an unidentified feature on the blue side, and 
is likely saturated. We derive a lower limit to its column density using the AOD technique. 

Figure 6 presents fits to some of the lines used in the analysis of the WD 1034+001 
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sightline. The line inside [] is of photospheric origin. Table 6 presents the column densities 
for D I, N I, 01, and Fe II along this sightline, obtained with the three different techniques. 

4.2.2. BD+393226 

Bluhm et al. (1999) report the detection of a weak interstellar absorption component 
at v & = -75 km s^ 1 (Component B in their Fig. 1), seen in 01 (A1039), Sill (A1526), and 
Fell (A2374). We see no evidence for such a component in the higher resolution and higher 
S/N FUSE data. Figure 7 displays the absorption profile of Fell 1144.9 as seen with FUSE. 
This transition is ~1.6 times stronger than the Fell A2374 transition in which Bluhm et al. 
(1999) claim to see a second component. The FUSE data shows only one component at v & 
= —25 km s -1 , consistent with their Component A at v Q = —25 km s _1 . Close inspection 
of the 01 A 971.7 transition (stronger than 01 A1039) in our data also shows only one 
absorption component at the same velocity of Fell A1144.9. There is however an absorption 
feature detected in the FUSE data, shifted by ~ —50 km s _1 from O I A1039 (consistent with 
absorption at v G = —75 km s^ 1 ). This feature is due to absorption by stellar NIII A1038.988 
and A1039.0, at the radial velocity of the star. It is possible that the lower S/N IUE data 
used by Bluhm et al. (1999) have lead them to mistakenly identify coincident noise features 
with extra absorption components of Fell and Sill. 

Although scattered flux is generally low in the FUSE channels (Moos et al. 2000; Sahnow 
et al. 2000), there is some apparent residual flux, <4% of the continuum, in the core of the 
HI Lyman lines (see Fig. 2). Apparent residual flux is also present, albeit at a smaller 
level, below the Lyman break. This is probably due to the fact that this target was observed 
in histogram mode. For this type of observation the FUSE pipeline takes into account a 
background model which is scaled only by the exposure time. For instance, WD 1034+001 
was observed in time-tag mode. In this case the background model used by the FUSE 
pipeline is scaled to match the observed counts in the unilluminated regions of the detector, 
hence doing a better job of removing the scattered flux (note how the cores of the HI lines 
in the spectrum of WD 1034+001 in Fig. 1 have no residual flux). To properly determine 
column densities it is important to account for this residual flux, which is approximately 
constant within each channel, but varies between channels. Consequently, we removed this 
residual flux from all the channels, before performing any analysis of the data (a similar 
procedure has been used in other D/H analyses, see Hebrard et al. 2002, for an example). 

The large separation between stellar and interstellar absorption (~ 255 km s" 1 Bluhm 
et al. 1999) along this sightline is high enough to remove the need for using a stellar model 
to account for the shape of the stellar continuum in the vicinity of the D I lines due to H I 
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and He II stellar absorption. Hence, no stellar model is used in this analysis. 

We use the AOD technique to determine column densities for DI, NI, 01, and Fell, 
and place lower limits on the column densities of CII, CII*, C III, Sill, PII, and Arl. 

With the COG method we fit single-component Gaussian COGs to the measured equiv- 
alent widths of DI, NI, 01, and Fell. No other atomic species detected along this sightline 
have enough transitions to perform COGs. From the four COGs we derive 6-values in the 
range 6 = 5-6 km s -1 . 

With the PF technique we fit a single component to unsaturated lines of D I, N I, O I, 
and Fell. H 2 is also included in the fit, but in a separate absorption component. The DI 
lines used in determining iV(D I) (see Table 3) with the PF method are not blended with 
H 2 or stellar absorption (to the best of our knowledge), hence iV(DI) is not affected by the 
H 2 column density. 

Weak absorption by O VI A1032 is detected at the same velocity as the one of the low 
ionization interstellar metals. OVI A1038 is blended with stellar absorption by OVP* and 
is not used in our analysis. We use AOD and PF to determine N(0 VI) along this sightline. 

Fits to some of the lines used in the analysis of this sightline are presented in Figure 8. 
Table 7 presents the column densities for D I, N I, O I, and Fe II along this sightline, obtained 
with the three different techniques. 

01 Column Density 

To determine N(OT) along the BD+39 3226 sightline we analyze the 919.9 (log (/A) = 
-0.79), 925.4 (log (/A) = -0.49), and 974 A (log (/A) = -1.82) 01 lines. We derive log 
N(01) = 16.33 ± 0.08 with the AOD method for the 919.9 A transition; a similar value is 
obtained for OI A925.4. We did not use the AOD technique on the weak A974 transition, 
or measured its equivalent width, because this line is blended on the red side with H 2 (J 
= 2). Figure 9 presents the curve-of-growth for 01 along this sightline. From the COG 
we derive log iV(OI) = 16.31 ± [j;^, also in good agreement with the AOD-derived iV(OI) 
values from the 919.9 and 925 A transitions (see Table 7). Figure 10 compares two sets of 
SiC IB and SiC 2A data, in the vicinity of the O I A974 transition (data shifted vertically for 
clarity). Data plotted at the top were calibrated with version 3 of CalFUSE (V3), data at 
the bottom were calibrated with version 2.4 (V2.4), which was used throughout this work. 
The discrepancy between V2.4 SiC IB and SiC 2A data around 01 A974 disappears when 
CalFUSE V3 is used. Hence, for the O I A974 transition, we use the V3 data (no other 
discrepancies between lines were found when the two datasets were compared). With the PF 
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technique we determine iV(OI) (as described in §4.2.2) by assuming that D I, NI, 01, and 
Fell trace the same absorption component. Hence, the 6-value of each species is constrained 
by the 6-values of the other species in the absorbing gas. Determining N(0 1) with O I A974 
only leads to log iV(OI) = 16.82 ± while a similar fit with only 01 A919.9 and A925.4 
leads to log AT(OI) = 16.51 ± 0.04. Fits using either A919.9 or A925.4 alone lead to JV(OI) 
consistent with N(0 1) quoted above when these two lines are fit together. 

The discrepancy between the N(0 1) values derived with or without the A974 O I line 
(> 0.30 dex in the log) could in principle indicate that there was a problem with saturation 
with some of the lines used in the analysis. However, the AOD determined N(0 1) for A919.9 
is in good agreement with the one determined with the same method from the A925 O I line, 
which is ~2 times stronger. This indicates that there is little or no unresolved saturation 
in these lines, and hence that N(0 1) derived from them is reliable. In addition, the COG- 
derived O I along this sightline (see Figure 9) is in good agreement with the AOD results 
quoted above. 

Overplotting log N(Ol) = 16.82 on the profiles of the A919.9 and A925.4 O I lines clearly 
indicates that this column density is overestimated, unless one assumes that some unknown 
effect causes simultaneously the 925.4 A and 919.9 A lines, in both channels, to be narrower 
and shallower than they should be. This seems unlikely. Note that the 6-values obtained in 
the fits of the three lines (A974 alone, or A919.9 together with A925.4, or all three transitions 
simultaneously) is similar. 

Three possibilities can be explored to try to explain the disagreement between N(0 1) 
derived from A974 and the one derived from the A925 and A919.9 OI lines: 1) there could 
be a problem with the /-value of the A974 transition, 2) there could be a problem with the 
/-values of A919.9 or A925, and 3) there could be a stellar line blended with the A974 OI 
transition and/or fixed pattern noise in one of the SiC channels (see below). A problem with 
the /-values of the A974 or A919.9 transitions seems unlikely as Hebrard & Moos (2003) 
have reported that no inconsistencies were found between A974 and A919.9 in the analysis 
of the BD+28 4211 sightline. In their OI analysis of the WD 2211-495 sightline Hebrard 
et al. (2002) performed tests in which N(0 1) derived independently from each individual O I 
transition (including amongst others the A919.9 and A925.4 transitions) was compared to 
N(0 1) derived from fitting all the O I lines simultaneously. They found that N(0 1) obtained 
from fits to the individual lines varied at most by 20% from the value obtained when all the 
lines were fit simultaneously and concluded that the /-values of the O I transitions used in 
their study were consistent with each other. The study by Hebrard et al. (2002) indicates 
then that there are no large discrepancies between the /-values of the A919.9 and A925.4 O I 
transitions (at most 20%). The discrepancy between iV(OI) derived from A974 and the one 
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derived from the A925 and A919.9 O I lines is much larger than 20%. Figure 10 compares the 
SiC IB and SiC 2A data for the BD+39 3226 sightline in the 01 A974 region. The positions 
of two neighboring H 2 lines are also marked. Even though O I is blended on the red side 
with H 2 ( J — 2), this blend should have little effect on the determination of N(0 1) since the 
two lines are resolved and N(H.2), J = 2, is well constrained by other weak unblended J = 
2 lines. The most likely possibility is then that part or all of the weak feature falling at the 
position of O I A974 along the BD+39 3226 sightline is of stellar origin. The spectra of the 
sdO BD+39 3226 contains many absorption lines, likely of stellar origin, that we are not able 
to identify, and so it is not surprising that one of them would fall at the expected position of 
the O I A974 line. We note that because stellar absorption is shifted by ~ —255 km s _1 from 
its rest frame wavelength, the possible blend of O I A974 with a stellar line does not affect 
the JV(OI) determination along other sightlines when the A974 transition is used. Finally, 
there is also the possibility, although unlikely, that N(0 1) derived from the A974 transition 
is correct and some unknown random effect causes iV(OI) derived from A919.9 and A925.4 
to agree even though these lines might suffer from saturation effects. Taking into account 
A r (0 1) obtained with the different methods and disregarding the value derived when the A 
974 transition is used, we adopt log A(OI) = 16.40 ± 0.10. 

Hebrard et al. (2005) and Friedman et al. (2005) have cautioned against not using the 
A974 O I line when deriving N(0 1) along high column density sightlines. In particular, 
Hebrard et al. (2005) has determined a new JV(OI) along the FeigellO sightline using the 
transition above and found it to be a factor of ~2 larger than the original N(0 1) reported 
by Friedman et al. (2002), which did not use the O I A974 line. There are however differences 
between our work and that of Hebrard et al. (2005) and Friedman et al. (2005) as neither 
of these authors used the apparent optical method in their analyses and in particular to 
compare N(Ol) derived from the A919.9 and A925.4 lines. As three of the authors of this 
paper also participated in the original analysis of the FeigellO sightline by Friedman et al. 
(2002) we have access to the same dataset used by Hebrard et al. (2005) in the new analysis. 
Using this dataset we performed AOD measurements of the A919.9, A925.4, and A974 OI 
lines and found that JV(OI) derived from A925.4 (log JV(OI) = 16.43 ± 0.03) is ~5cr smaller 
than iV(OI) derived from the A919.9 line (log JV(OI) = 16.58 ± 0.03). From the A974 line 
we derive log N(OT) = 16.98 ± 0.08 (la), in perfect agreement with the new log N(OT) = 
17.06 ± 0.15 (2a) quoted by Hebrard et al. (2005). The different results we obtained from 
the A919.9 and A925.4 lines for the Feige 110 sightline are a strong indication of saturation of 
one or both of these two lines, and unless one compares the AOD derived iV(OI) from A919.9 
with AT(OI) from a weaker line one can not decide whether A919.9 suffers from saturation or 
not. As mentioned before, in our case the N(Ol) results from A919.9 and A925.4 lines are 
in excellent agreement, which gives us confidence that our results are not biased due to not 
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using the A974 transition. 

4.2.3. TD 132709 

Along this sightline the separation between stellar and interstellar absorption is ~13 km 
s -1 . We use the stellar model computed to determine AT(HI) (see §4.3 below) to ascertain 
which ISM lines might suffer from stellar blendings. The stellar model indicates that the 
profiles of the D I lines used in this analysis (see Table 3) are not affected by the shape of the 
stellar continuum in the vicinity of these lines, hence no stellar model is needed to determine 
the column densities discussed below. In addition, the small molecular hydrogen content 
along this sightline (log iV(H 2 ) = 14.48 ± on, see §4.1.3) implies that the deuterium column 
density derived below is not sensitive to iV(H 2 ). 

We use the AOD technique on weak lines to determine the column densities of D I, 
NI, 01, and Fell, which are then compared to iV derived with the other methods. CII, 
CTI*, Sill, PII, and Arl only have saturated transitions in the FUSE bandpass. We use 
the AOD method to determine lower limits to their column densities. We use the COG 
technique to determine the column densities of DI, NI, 01, and Fell. The b- values derived 
from the four COGs fall in the range 3-6 km s _1 . With PF we fit a single absorption 
component to unsaturated lines of DI, NI, 01, and Fell. H 2 is also included in the fit, but 
in a separate absorption component. In order to provide a smooth continuum for the DI 
lines we also fit H I as an independent component. This fit, however, is not used to determine 
iV(HT). Absorption by O VI (A1032, 1037) is also detected along this sightline but, the small 
separation between stellar and interstellar absorption along this sightline (~ 13 km s _1 ) does 
not allow us to determine if this absorption is of stellar or interstellar absorption origin (or 
both), as stellar and interstellar absorptions are blended. 

The curve of growth for O I along this sightline is presented in Figure 11. We derive log 
N(0 1) = 16.45 ± QQ3. Figure 12 presents fits to some of the lines used in the analysis of the 
TD 132709 sightline. Table 8 presents the column densities for DI, NI, 01, and Fell along 
this sightline, obtained with the three different techniques. 
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4.3. HI Column Densities 

4.3.1. WD 1034+ 001 

To determine iV(HI) along this sightline we use data obtained by the GHRS onboard 
HST. GHRS observations were obtained in June of 1992, with the G160M grating and 
the LSA aperture (R ~ 20,000) for a total exposure time of 288 s (rootname: z0ye0c08t). 
Unfortunately these observations cover only the wavelength range 1185 - 1221 A, which 
does not include the complete red wing of Lyct. However, the existing data still allow us to 
determine iV(HI). 

The calibrated data were retrieved from the Multimission Archive at the Space Telescope 
Science Institute and no further processing was applied. The error-weighted data, obtained 
at four different FP-split positions, were combined to increase the signal to noise ratio of 
the final dataset. We consider the influence on the Lya profile of several stellar models with 
different effective temperatures and gravities. These are described in detail below. 

Stellar Model 

We computed stellar atmosphere models using the atmospheric parameters determined by 
Werner et al. (1995), i.e., T eff = 100,000 ± ^ooo K > lo g 9 = 7 - 5 ± °- 3 > and lo g ( He / H ) = 
3.0. NLTE H+He models were computed with the stellar atmosphere codes TLUSTY and 
SYNSPEC (see, e.g., Hubeny & Lanz 1995). The line profile of the He IIA1215 line uses 
an approximate Stark broadening treatment developed by Hubeny et al. (1994), while the 
hydrogen line profile uses the Vidal et al. (1973) theory (Lemke 1997). We define as the 
best fit model the stellar model with T e g = 100,000 K, log g = 7.5, and log (He/H) = 3.0. 
The 1 a uncertainties in the stellar parameters determined by Werner et al. (1995) are used 
to compute stellar models that allow us to determine the systematic uncertainties in iV(H I) 
associated with the uncertainties in the stellar models. Determining the H I column density 
with different stellar models places more credible error bars on iV(HI) than using a single 
stellar model. The statistical uncertainties are derived by using the best fit stellar model. 
The systematic uncertainties associated with the stellar models are estimated by considering 
the most extreme stellar models, i.e., those that yield the strongest and weakest stellar Lya. 
Figure 13 illustrates such stellar profiles. The figure shows that the best fit stellar model 
(T e ff = 100,000 K, logg = 7.5; solid line) produces a Lja line profile that is between the 
strongest Lya line (T eS = 90,000 K, \ogg = 7.8; dash-dotted line) and the weakest one 
(T cff = 115,000 K, logs = 7.2). 
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Lya Profile Fitting 

We use the separation between stellar and interstellar absorptions derived by Sion et al. 
(1985) and Werner et al. (1995) to align the data and the stellar models. The different 
models are then scaled by different factors in order to match the spectrum flux over the 
region 1190 - 1195 A. The spectrum is then normalized by the different models prior to 
fitting. 

To determine the H I column density we fit three absorption components with indepen- 
dent velocities, one each for HI, DI, and H 2 (J = - 3). N(Dl) and iV(H 2 ) are fixed at 
the adopted values (see Tables 5 and 9) during the fitting procedure. The 6-value of the 
component containing DI is fixed at b = 5.8 km s -1 , derived from the DI COG. In order 
to constrain 6 HI we fit also the higher order HI Lyman lines 917 A, 918 A, 919 A, 920 A, 
and 926 A (923 A is blended with stellar NIV). The Si III line (A1206) is also included in 
the Lya fit, in a fourth independent absorption component, so we can use a large continuum 
region in the blue side of the Lya wing. The normalized data is fitted in the range 1201.15 
- 1221.28 A; the continuum is fixed at 1.0 during the fitting procedure. 

The statistical uncertainties associated with iV(HI) are determined by fitting the data 
normalized by the best fit model, following the procedure outlined in Oliveira et al. (2003). To 
determine the uncertainties associated with the continuum placement due to the uncertainties 
in the stellar model we perform two additional fits, in which the continuum level is multiplied 
by factors of 1.07 and 0.93 during the fitting procedure (corresponding to a 7% change in 
the continuum level or to a by-eye' of 2 a). As mentioned above, the effect of systematic 
uncertainties associated with the stellar models are determined by using the extreme models 
discussed above to measure the highest and lowest interstellar iV(HI). These models are 
illustrated in Figure 13. Figure 14 presents the fit to the Lya region along this sightline, 
using the best fit stellar model to normalize the data. 

We take a conservative approach to combine the different uncertainties by adding them, 
rather than adding them in quadrature. Taking into account the different uncertainties in 
the manner discussed above we determine log iV(HI) = 20.07 ± 0.07 (lcr). 

4.3.2. TD 132709 

To determine iV(H I) along this sightline we use data obtained with /[/Eand FUSE. High 
dispersion, large aperture IUE observations of the Lya region were obtained in March 1981 
(swpl3459), March 1994 (swp50226, swp50227), and December 1993 (swp49676, swp49677), 
for a total exposure time of 53940 s. The calibrated NEWSIPS MXHI data were retrieved 
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from the Multimission archive at the Space Telescope Science Institute. The data were not 
processed further. 

The different observations of Lya were coadded (weighted by their uncertainties) to 
increase the signal-to- noise ratio of the final dataset. A set of stellar models, described 
below, were computed to take into account the stellar absorption in the vicinity of Lya. 
Similarly to WD 1034+001, we consider the influence of several stellar models with different 
temperatures and gravities, on the Lya profile. 

Stellar Model 

We computed NLTE H+He+C+N models using the atmospheric parameters obtained by 
Dreizler (1993) from fitting the optical spectrum of TD1 32709. These authors determined 
T eff = 46,500 ± 1,000 K, log g = 5.55 ± 0.15, and log(He/H) = 2. We define as the best fit 
model the model corresponding to T eff = 46,500 K, log g = 5.55, and log(He/H) = 2, and 
use the uncertainties in the stellar parameters to compute stellar models which are used to 
determine the uncertainties in iV(HI) associated with the stellar models. We considered also 
the influence on the adopted He/H by computing another model with T e g = 45,000 K, logg 
= 5.7 and log(He/H) = 1. This produced iV(HI) within the uncertainties of the adopted 
value. 

Figure 15 presents the best fit stellar model (solid line), corresponding to T c r = 46,500 
K and logg = 5.55, together with the models that produce the strongest (T e g = 45,500 K 
and logg = 5.7; dashed line) and weakest (T e ff — 47,500 K and logg = 5.4; dash-dotted 
line) stellar absorption. These models are used to determine the systematic uncertainties in 
iV(HI) associated with the stellar models. 

Lya Profile Fitting 

We use the stellar CIV doublet at A1230 to align the stellar models with the data. The 
stellar models are then scaled by different factors to match the spectrum flux in the 1200 
A and 1225 A regions. Prior to fitting we normalize the data by the different stellar models. 

To determine N(H I) we fit three absorption components, one each for H I, D I + O I, and 
H 2 (J = - 2). The column densities of all the species but H I, are fixed at the adopted values 
(see Tables 9 and 5) during the fitting procedure. The 6-value of the component containing 
DI and OI is fixed at b = 6.0 km s _1 , derived from the DI and OI COG analyses. Higher 
order Ly lines (916 - 919 A) are also included in the fit, to constrain &hi- These Ly lines 
were chosen because our models indicate that the interstellar absorption spectrum is not 
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affected by stellar absorption (either from H or He) in this region. Geocoronal emission 
fills part of the core of the Lyct line and is taken into account during the fitting procedure 
by introducing a background in the fit. We determine the background level using the blue 
side of the Lyct core; the airglow emission visible in Figure 15 is not allowed to influence 
the profile fitting procedure. We estimate the errors associated with this background by 
performing multiple fits, using the best fit model, where the background is varied between 
6 and 12% of the continuum level. Similarly, errors associated with the continuum level 
are determined by varying the continuum level, using the best fit model. Figure 16 (top 
panel) presents the fit to the Lya interstellar absorption when the best fit stellar model 
is used to normalize the data and with the zero-flux level defined by the blue side of the 
Lja wing, yielding log iV(HI) = 20.03. A similar fit, but with the zero-flux level defined by 
the red side of the Lya wing is shown in the bottom panel of this figure. This fit yields log 
iV(HI) = 20.12. The small emission line on the red side of the airglow emission is likely an 
artifact of the IUE data similar to the ones studied by Crenshaw et al. (1996). Taking into 
account this feature when defining the zero-flux level leads to JV(HI) in the range quoted 
above (20.3-20.12). The statistical uncertainties associated with iV(HI) are determined by 
fitting the data normalized by the best fit model, following the procedure outlined in Oliveira 
et al. (2003). Taking all these uncertainties into account, by adding them as in the case of 
WD 1034+001, we determine log iV(HI) = 20.03 ± 0.10. 



5. RESULTS 

Table 10 presents several ratios, using the column densities summarized in Table 9. 
The fraction of H 2 along these sightlines, / Ha = 2xiV(H 2 )/(2x iV(H 2 ) + iV(HI)), is log / Ha 
= -4.05 ± 0.58, -4.13 ± 0.46, and -5.25 ± 0.80 for WD 1034+001, BD+39 3226, and 
TD1 32709, respectively. Hence, for the three sightlines, JV(H) = AT(HI) + 2xA r (H 2 ) « 
iV(HI). 



5.1. D/H Ratios 

We derive (D/H)xl0 5 = 2.14 ± 1.17 ± ° H, and 1.86 ± for WD 1034+001, 
BD+39 3226, and TD1 32709, respectively. Only the D/H ratio toward TD1 32709 is con- 
sistent, at the la level, with the LB value (D/H = (1.56 ± 0.04) xl0~ 5 ) derived by Wood 
et al. (2004). The value of D/H toward WD 1034+001 is one of the highest D/H ratios 
measured in the nearby interstellar medium. The other high values are (1.91 ± q 2 4)x10~ 5 
(PG 0038+199, Williger et al. 2005), (2.14 ± 0.41) xl0~ 5 (FeigellO, Friedman et al. 2002), 
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(2.18 ± ° f g )xl0- 5 ( 7 2 Vel, Sonneborn et al. 2000), (2.24 ± x lO" 5 (a Cm, York & 
Rogerson 1976) and (2.24 ± ° f e )xl0' 5 (LSE44 Friedman et al. 2005). 

The D/H values reported above are averages over the entire sightline; however they 
are dominated by the gas beyond the Local Bubble. WD 1034+001 is the shortest sightline 
analyzed in this work, d = 155 ± H pc. Even though the uncertainties on its distance 
put it close to the Local Bubble boundary (~100 pc, Sfeir et al. 1999, however, in some 
directions the boundary can be at 65 or 250 pc), the bulk of the gas probed by this sightline 
is outside the Local Bubble (log iV(H) < 19.2). The Local Bubble contribution to the 
WD 1034+001 D and H column densities can be estimated by assuming logiV(H) = 19.2 and 
using D/H = 1.56xl0 -5 derived by Wood et al. (2004) for the Local Bubble to determine 
logiV(D) L B = 14.39. Subtracting the LB contribution yields then D/H = 2.23xl0~ 5 for the 
gas outside the LB, an insignificant change from the D/H ratio quoted in Table 10 for this 
sightline and consistent with the quoted uncertainties (D/HxlO 5 = 2.14 ± o;||). For the 
BD+39 3226 and TD1 32709 subtracting the LB contribution to the sightline D/H ratios 
in the manner described above leads to D/HxlO 5 = 1.12 and 1.92, respectively, consistent 
within the uncertainties with the ratios reported in Table 10 for these two sightlines. 



5.2. O/H Ratios 

We derive (O/H)xl0 4 = 3.39 ± 2.09 ± ° f 8 , and 2.45 ± for WD 1034+001, 
BD+39 3226, and TD1 32709, respectively. The first and last ratios are consistent, within 
the quoted uncertainties, with O/H = (3.43 ± 0.15) x 10 -4 derived for 13 lines of sight probing 
gas within 1500 pc (most within 500 pc), with 20.18 < log JV(HI) < 21.28 (Meyer et al. 1998, 
updated /-value). O/H along the BD+39 3226 sightline is »1.86<t away from the Meyer et al. 
(1998) ratio. Taking into account the Local Bubble contribution to the O/H ratio along the 
three sightlines using logiV(H) = 19.2 and the O/H ratio derived by Oliveira et al. (2005) 
for the Local Bubble, O/Hxl0 4 = 3.45 ± 0.19, leads to O/H ratios for the gas beyond the 
LB consistent within the uncertainties with the values quoted in Table 10. 



5.3. D/O Ratios 

We derive (D/O)xl0 2 = 6.31 ± \f 8 , 5.62 ± \*\, and 7.59 ± 2 ;^, for WD 1034+001, 
BD+39 3226, and TD1 32709, respectively. Hebrard & Moos (2003) derived D/O = (3.84 ± 
0.16)xl0~ 2 for the Local Bubble from measurements along 14 sightlines. All the D/O ratios 
reported here are inconsistent, at the la level, with the Local Bubble value, and are some 
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of the highest D/O ratios measured in the ISM. The implications of this will be discussed 
below. 

5.4. N/H, O/N, and D/N Ratios 

The N/H ratios computed for the three sightlines, (N/H)xl0 5 = 7.76 ± 2 2 f , 5.89 ± 
and 8.91 ± || (for WD 1034+001, BD+39 3226, and TD1 32709, respectively), are all con- 
sistent, within the uncertainties, with N/H = (7.5 ± 0.4)xl0~ 5 derived by Meyer et al. 
(1997). 

We derive O/N = 4.37 ± \ f 8 , 3.55 ± \ f 3 , and 2.75 ± J;gJ. These are consistent, at the 
la level (1.34a for TD1 32709), with O/N = 4.1 ± 0.3 derived using the values of O/H and 
N/H determined by Meyer et al. (1998) and Meyer et al. (1997). 

The D/N ratios for the three sightlines are (2.75 ± ^xlO -1 , (2.00 ± g;^)xl0 _1 , and 
(2.09 ± ^xHT 1 (for WD 1034+001, BD+39 3226, and TD1 32709, respectively). 

5.5. N from FUSE vs. N from ORFEUS for the BD+39 3226 sightline 

Table 11 displays, side by side, the column densities of several species derived in this work 
from FUSE data and the column densities derived by Bluhm et al. (1999) using ORFEUS 
and IUE data. The column densities for DI, OI, Fell, and H 2 (J = 0, 1, 3) are consistent 
within the quoted uncertainties. For these species, using the FUSE data to determine their 
column densities leads in all cases to a significant improvement of the uncertainties. This is 
particularly important for D I, N I, and O I, for which we want to determine ratios relative 
to HI as accurate as possible. The value of iV(NI) determined by Bluhm et al. (1999) 
is substantially smaller than the one determined in this work; log iV(NI) = 14.75 ± 0.25 
vs. 15.85 ± 0.10. We suspect this is because the weakest NI line used by Bluhm et al. 
(1999) to determine iV(NI) (A964.6256) suffers from considerable saturation, leading to an 
underestimation of the value iV(NI). 
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6. DISCUSSION 

6.1. WD 1034+001 

The WD 1034+001 sightline has one of the highest D/H ratios reported in the literature, 
D/H = (2.14 ± 0;|)xl0- 5 , 1.29a above D/H LB derived by Wood et al. (2004). However, 
its O/H and N/H ratios are consistent at the la level with previous measurements (Meyer 
et al. 1998, 1997). A D/H ratio consistent with the LB value (~1.56xl0 -5 ), would require 
logA(HI) = 20.21. In such case, the O/H and N/H ratios (2.47 ± g;i§)xl0 _4 and (5.66 
± i"0q)x1O -5 , respectively) would still be consistent with the values derived by Meyer et al. 
(1998) and Meyer et al. (1997), while not affecting the high D/O ratio (the same percentage 
uncertainties were used for the required iV(HI) as the ones derived in this work). In this 
scenario ~ 34% of the HI would have been missed in our analysis. It seems unlikely that 
iV(HI) along this sightline could be underestimated by such a large factor. On the other 
hand, a value of logiV(D I) = 15.26 (38% smaller than our derived value) together with our 
derived iV(HI) would also lead to D/H consistent with the LB value, without affecting the 
O/H and N/H ratios. In this scenario, using the same uncertainties for the required iV(DI) 
as the ones derived in this work, we would find D/O = (4.57 ± J^g)xl0 -2 . This is lower 
than our derived value of (6.31 ± J'^xlO -2 , and consistent with the Local Bubble value 
derived by Hebrard & Moos (2003). We derived iV(DI) along this sightline using three 
different methods, AOD, COG, and PF (see Table 3), and it is unlikely that AT(DI) has 
been overestimated. 

Linsky et al. (2005, in prep) looked at correlations between the D/H ratio and other 
parameters, such as the depletions of Fe and Si, and the kinetic temperature of H 2 (T i), in 
order to try to understand what is the total deuterium abundance in the local Galactic disk. 
They found that sightlines with higher Fe and Si depletions correspond in general to lower 
D/H ratios and that sightlines with high T 01 have typically higher D/H ratios. 

The Fe/H ratio for this sightline, Fe/H = (1.07 ± gp)xl0~ 6 (corresponding to log 
(Fell/H I) = —5.97), is lower than would be expected for this D/H ratio according to Figure 
3 of Linsky et al. (2005, in prep). Using the solar abundance of Fe from Asplund et al. 
(2004), [Fe/H] = -4.55, the depletion of Fe along this sightline is D(Fe) = -1.42. This is 
an unusual sightline, for which non-stellar absorption by C IV, N V, O VI, Si IV, S III, and 
S IV, has been detected (Sion et al. 1985; Werner et al. 1995, and this work), in addition to 
the large elliptic halos seen in emission (see §2.1). 

The presence of highly ionized species, the high D/H and lower than expected Fe/H, 
and the typical O/H and N/H ratios are consistent with a scenario in which D is released 
from the grains, but not much Fe. Perhaps the high radiation field of WD 1034+001 plays 
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a special role for this sightline. A similar mechanism has been proposed by Linsky et al. 
(2005, in prep) to explain the high D/H and low Fe/H ratios along the sightline to 7 2 Vel. 

6.2. BD+39 3226 

For the BD+39 3226 sightline the D/H ratio is only slightly inconsistent with the LB 
value (1.15a below). For this sightline, D/H and Fe/H ((1.17 ± o:i) xlo ~ 6 ) follows the 
correlation derived by Linsky et al. (2005, in prep) for these quantities. Even though 
interstellar O VI is detected along this sightline, we found no evidence of other highly ionized 
species (such as C IV, N V, Si IV, S III, or S IV) in either the FUSE or IUE data for this star. 
In this case the separation between stellar and interstellar absorption is high enough (~255 
km s _1 ) to allow for a definitive detection, should an interstellar feature be present. It 
is possible that along this sightline OVI is formed in evaporative interfaces between cool 
clouds and the hot and diffuse gas in the Local Bubble. This mechanism has been invoked 
by Oegerle et al. (2005) to explain the presence of O VI along the sightlines to several white 
dwarfs inside the Local Bubble. D/O along this sightline, (5.62 ± \-%\) x 10 -2 , is not consistent 
at the 1 a level with the D/O ratio derived for the LB by Hebrard & Moos (2003), D/O = 
(3.84 ± 0.16) xl0~ 2 . This is however, in part due to the low O/H ratio derived here for this 
sightline, O/H = (2.09 ± q^xIO -4 , which is 1.86(7 below the O/H ratio derived by Meyer 
et al. (1998). This O/H ratio is lower than what has been derived in similar studies (see 
§7.1) and might be indicative of variations of the O/H ratio. Using the O/H ratio derived 
by the authors above yields D/O ~ 3.4xl0~ 2 , in closer agreement with the D/O LB value 
reported above. 

6.3. TD1 32709 

The TD1 32709 sightline presents a D/H ratio consistent at the la level with the LB, 
while D/O is inconsistent within the quoted uncertainties at the la level. O/H is lower 
than the value derived by Meyer et al. (1998) by 1.09a. If we assume that the D/H and 
D/O ratios are high due to iV(DI) being overestimated then we would need logiV(DI) = 
15.22 to bring D/H closer to the LB value. Using the same uncertainties for the required 
DI as the ones derived in this work (0.05 dex), the D/O ratio would still be high at (6.31 
± x 10 -2 , even with this reduced value for iV(DI). It is unlikely that iV(HI) has been 
over or underestimated, as a lower iV(HI) would lead to an even larger D/H ratio, while 
a larger iV(HI) would lead to an even smaller O/H ratio. Increasing iV(OI) so that O/H 
agrees with the Meyer et al. (1998) value still leads to a high D/O ratio, (5.43 ± J;^) xl0~ 2 
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(assuming that the increased N(0 1) has the same uncertainties as the one derived in this 
work). The only possibility to reconcile the D/H and D/O ratios along this sightline with the 
LB values would be for D and O to be simultaneously over and underestimated by 1.6 a and 
1.5 a, respectively. However this is unlikely given the number of transitions and methods 
used to determine iV(DI) and iV(OI). 

For this sightline the Fe/H ratio ((0.83 ± 24) x 10~ 6 ) is lower than would be expected 
for the D/H ratio according to Figure 3 of Linsky et al. (2005, in prep). We searched the 
IUE data for this star for absorption by highly ionized species such as C IV, N V, and Si IV. 
However because the separation between stellar and interstellar absorption is low (~13 km 
s _1 ) we cannot determine if the detected absorptions are of stellar or interstellar origin (or 
both), similarly to what we had concluded for the OVI absorption detected in the FUSE 
data (see §4.2.3). 

6.4. Implications for models that explain D/H variations 

The high D/H ratios derived for the WD 1034+001 and TD1 32709 sightlines increase 
the sample of sightlines for which a high D/H ratio has been found. These measurements 
strengthen the argument that these high ratios can not be due to some systematic problem 
associated with the different analyses (see also discussion below) or that they are peculiar 
cases. In particular, it has been argued that the high D/H ratios found in the ISM could be 
due to systematic effects that affect the jV(H I) measurements (Hebrard et al. 2005). Because 
the D/O ratios derived for these two sightlines are also high and independent of iV(HI), the 
high D/H ratios are likely due to a high abundance of D and not to a systematic problem 
with the iV(HI) determination (see also discussion in §7.1). In addition, the high D/O 
ratios, together with O/H ratios consistent with the Meyer et al. (1998) results, contradict 
the argument by Hebrard (2006) against a high present-epoch D/H ratio, on the basis that 
no correspondingly high D/O ratios have been measured. It is also hard to explain these high 
D/H ratios on the basis of localized infall of metal poor gas, as one would expect the O/H 
and N/H ratios for these sightlines to be affected as well. In addition, several studies of the 
O/H ratio (Meyer et al. 1998; Andre et al. 2003; Cartledge et al. 2004; Oliveira et al. 2005) 
have shown that this ratio is constant across a large range of N(0 1), distance to targets, and 
average sightline gas density. However, note that in a special case where the infalling gas had 
a high D/H ratio and metallicity comparable to the local ISM value of Meyer et al. (1998), 
such a mechanism could not be ruled out. Variable astration in which the amount of D burned 
in the interior of stars and O produced, varies from sightline to sightline, seems also to be 
an unlikely mechanism for explaining the high D/H ratios. In the case of WD 1034+001 and 
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TD1 32709 less deuterium would have to have been destroyed by astration and consequently 
less O would have been produced by supernovae. This scenario seems implausible because 
it implies that astration rates would have to be variable on short distance scales in order to 
explain the nearby sightlines with LB-like D/H ratios and also because it is not supported 
by our derived O/H and N/H ratios. In addition, in both the infall and variable astration 
scenarios, these processes would have to occur in a time smaller than the typical mixing 
time-scale of 350 Myr (de Avillez & Mac Low 2002). Our results support the idea that 
some sightlines in the Milky Way ISM have high D/H ratios. Linsky et al. (2005, in prep) 
have used the existence of these sightlines together with the deuterium dust-depletion model 
of Draine (2004) to argue for a high present-epoch abundance of deuterium. Linsky et al. 
(2005, in prep) estimate the total abundance of D (in gas + dust) in the local disk of the 
Galaxy, D/H > (2.19 ± 0.27) xl0~ 5 , by taking the weighted average of the ratios for the five 
sightlines with the highest values. Our new measurements for WD 1034+001 and TD1 32709 
are in good agreement with this estimate. Below we compare our new measurements with 
ratios from the literature. 

7. COMPARISON WITH PREVIOUS MEASUREMENTS 

7.1. D/H, O/H, and D/O vs. iV(H) 

Figure 17 presents D/H, O/H and D/O ratios (lcr uncertainties, top, middle, and bottom 
panels, respectively) as a function of the total hydrogen column density, iV(H), along the 
sightline. The data used in this figure are presented in Table 12. Sightlines for which all the 
ratios are available are marked by asterisks (literature values) and filled circles (this work). 
The data in this plot uses the ratios derived in this work together with the compilation 
by Wood et al. (2004) and the values of Friedman et al. (2005) and Hebrard et al. (2005). 
Sightlines for which no N(0) is available are marked with open squares in the top panel 
(uncertainties displayed as dotted lines for clarity of the plot). Sightlines displayed by open 
squares with dashed error bars in the bottom panel (corresponding to D/O) do not have 
iV(H) measurements available. For these sightlines we estimate iV(H) using iV(OI) and 
the O/H ratio of 3.43xl0" 4 derived by Meyer et al. (1998). Also plotted are the Local 
Bubble values of D/H = (1.56 ± 0.04) xl0~ 5 (top panel, Wood et al. 2004) and D/O = 
(3.84 ± 0.16) xl0~ 2 (bottom panel, Hebrard & Moos 2003). Dashed vertical lines mark the 
approximate position of the Local Bubble (log N(R) < 19.2), and log N(H) = 20.7, where 
a new low and relatively constant D/H ratio seems to emerge. The value of O/H = (3.43 ± 
0.15) x 10 -4 derived by Meyer et al. (1998) for the local ISM is also marked (solid and dashed 
lines in the middle panel of Figure 17). 
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For the range of hydrogen column densities 19.2 < log iV(H) < 20.7, assuming a constant 
D/H ratio yields the weighted mean (D/H)xl0 5 = 1.18 ± 0.05 (1 a on the mean) and xl 
= 5.9 for 16 degrees of freedom (where xl is the total x 2 divided by the number of degrees 
of freedom, v). Because some of the ratios in Table 12 have asymmetric uncertainties we 
use the average of the uncertainties on each ratio to compute the weighted mean and the 
uncertainty in the weighted mean for D/H and for the other ratios discussed below. One can 
consider two possibilities to explain the high value of xl- The scatter is real and variations 
of the D/H ratio do exist, or the uncertainties of the individual measurements have been 
underestimated. This last possibility can be explored by assuming that the scatter in D/H 
for 19.2 < logiV(H) < 20.7 is due to some random systematic error in the determination of 
D/H, unaccounted for in the quoted uncertainties. One can then estimate by how much the 
uncertainties would have to be increased to obtain xl — 1- Assuming further that this error 
can be expressed as a fraction of D/H, and that it is added in quadrature to the quoted 
uncertainty 0d/h ; <^ new = (c^/h + f 2 x(D/H) 2 ) ' 5 , we derive / = 0.40. This implies that the 
unknown systematic errors would have to be 40% of D/H. Such a high uncertainty seems 
unlikely. If we also add the ratios for logiV(H) > 20.7 to the dataset, we derive the weighted 
mean D/H = (1.10 ± 0.04) xl0~ 5 with xl = 5.1 for 21 degrees of freedom, and calculate / = 
0.34, still a high and unlikely value. These results show that it is unlikely that the scatter in 
D/H is due to unknown systematic uncertainties associated with the different measurements, 
and that it is real. 

As mentioned above, the apparent lack of scatter in the D/O ratios, has been used as 
an argument against the existence of sightlines with high D/H ratios (Hebrard et al. 2005). 
One has however to be careful in analyzing these ratios, as they depend on both the N(D) 
and N(0) values. For instance, there are different ways of producing a low D/O. One can 
have a low D/H ratio, and normal O/H, or a normal D/H and a high O/H (where normal 
refers to the LB value for D/H and to the Meyer et al. 1998, value for O/H). This last way 
must not be common as several studies have concluded that O/H is constant throughout a 
large range of iV(H) (see for example Andre et al. 2003; Cartledge et al. 2004, and Oliveira 
et. al. 2005). However, as shown in Figure 17, four out of the five sightlines with log iV(H) > 
20.7, also have O/H ratios inconsistent with the Meyer et al. (1998) value (solid and dashed 
horizontal lines in the middle plot). For the mean of these five points we find O/H = (5.33 
± 0.36) xlO -4 . However xl = 5-3 for four degrees of freedom, indicating that a single mean 
may not be appropriate. In other words, for four of these targets D/O is lower than it would 
have been, if these sightlines had normal O/H ratios. Because of the high average O/H 
value and considering that the number of measurements at large A(H) (>20.7 in the log) 
is statistically small, we think it is necessary when estimating the deuterium abundance at 
high column densities (which according to the suggestion of Hebrard & Moos 2003, may be 
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the present epoch Milky Way abundance) to consider both D/H and D/0 measurements. In 
addition, the use of the Meyer et al. (1998) value for 0/H to convert the average D/0 in this 
region to the average D/H yields an estimate that is low when compared to the direct D/H 
measurements in this iV(H) region, not unlike what happens when we do the same exercise 
for the Local Bubble. 

The five points with high H column density presented in Figure 17 are near the limit 
of deuterium column densities that can be measured with an instrument with FUSE-like 
properties. The deuterium lines used in these analysis (Hoopes et al. 2003; Wood et al. 
2004; Hebrard et al. 2005) would have been saturated if the D/H ratios along these sightlines 
were similar or larger than the LB D/H value (~1.56xl0~ 5 ). Hence, there is a bias against 
measuring high D/H ratios at these large H column densities. This could explain in part 
why only a few high D/H ratios have been reported in the literature. 

Additionally, to try to understand the behavior of the D/O ratio in the regime 19.2 < 
logiV(H) < 20.7 and the apparent lack of scatter in D/O compared to D/H, we can consider 
the O/H ratios for the sightlines for which D/H > 1.9xl0~ 5 in the top panel of Figure 17 (a 
Cru, 7 2 Vel, WD 1034+001, FeigellO, PG 0038+199 and LSE44; as a function of increasing 
iV(H)). For these sightlines with high D/H, one would expect also a high D/O, assuming 
that O/H is similar to the Meyer et al. (1998) value, hence producing more scatter in the 
D/O measurements. However, two of these sightlines do not have O/H measurements (a Cru 
and 7 2 Vel), while FeigellO, PG 0038+199, and LSE44 have abnormally high O/H ratios 
(see Table 12). The remaining sightline (WD 1034+001) has O/H consistent with the Meyer 
et al. (1998) value. Hence, the apparent lack of scatter in the D/O ratio, which has been 
used as an argument against a high present-epoch D/H ratio, is the result of small number 
statistics, and cannot be used to draw definitive conclusions about the behavior of the D/O 
ratio in this regime of iV(H). 



7.2. D/O vs. Fe/O 

Linsky et al. (2005, in prep) have looked at the correlation between the deuterium 
abundance and the depletion of Fe and found that in sightlines with low D/H ratios Fe 
is typically more depleted than in sightlines with higher D/H ratios. A different way of 
considering the relationship between the deuterium and iron abundances is to look at D/O 
versus Fe/O. By taking the ratio of these abundances relatively to oxygen rather than to 
hydrogen one can avoid systematic problems that might affect the determination of iV(H I) 
(as could be the case for the FeigellO sightline, as argued by Hebrard et al. 2005). Figure 
18 presents D/O as a function of Fe/O. The solid line represents the fit to all data points 
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in this plot, yielding D/O = (2.04 ± 0.15)xl0~ 2 + (1.56 ± 0.30)x 10- 2 xFe/O, with X 2 = 
76 for 18 degrees of freedom (only uncertainties in D/O taken into account). This fit is 
mostly influenced by the sightline with the largest Fe/O ratio and small D/O uncertainties, 
G191— B2B, represented by a square in the plot. Removing this sightline from the fit yields 
D/O = (0.99 ± 0.23)xl0- 2 + (5.64 ± 0.73)xl0~ 2 xD/O, with X 2 = 39 for 17 degrees of 
freedom. The other three sightlines represented by squares also have Fe/O ratios which 
are inconsistent with the trend displayed by the majority of the data plotted in Figure 18. 
Excluding all the sightlines displayed by squares in the plot from our fit increases the slope 
slightly more yielding D/O = (0.47 ± 0.28) xlO" 2 + (7.85 ± 0.99) xl0~ 2 x Fe/O, with X 2 
= 16 for 14 degrees of freedom. The probability of X 2 being larger than 16 for 14 degrees 
of freedom is ~30%. Thus, the trend of decreasing Fe/H with decreasing D/H observed by 
Linsky et al. (2005, in prep) is also observed when one considers D/O versus Fe/O. We note 
however that there is no particular reason for not taking into account the points displayed 
in Figure 18 by squares or that it is theoretically expected that N(D) depends linearly on 
AT(Fe). 



7.3. D/H and (n H ) 

Another way of trying to understand the behavior of D/H in the ISM, besides looking 
for correlations with the abundances of other species (e.g., Prochaska et al. 2005, Linsky et 
al. 2005, in prep), is to study D/H as a function of the average sightline gas volume density, 
(n H ) = N(R)/d. 

Elemental abundances measured in the interstellar gas in the solar vicinity are generally 
depleted with respect to the solar values (Morton 1974; Jenkins et al. 1986). This depletion 
effect has been commonly interpreted as the result of a fraction of the elements being locked 
up in dust grains (Jenkins 2004). Furthermore, several studies have shown that the strength 
of the depletions increases as the average gas density increases (see e. g., Savage & Bohlin 
1979; Harris et al. 1984; Spitzer 1985; Jenkins et al. 1986; Jenkins 1987; Crinklaw et al. 1994). 
The idea that D might be depleted in grains was first proposed by Jura (1982) and updated 
recently by Draine (2004) (see also Draine 2003). Harris et al. (1984) looked for a correla- 
tion between the deuterium abundance and the sightline density. Using Copernicus data 
for 14 sightlines and a limited range of densities, they found no evidence of significant deple- 
tion of interstellar deuterium. However, since Copernicus, several space-born observatories 
(I MAPS, HST, FU SE) have been used to measure the interstellar abundance of deuterium 
and there are now D/H measurements for more than forty sightlines. It is then appropriate 
to revisit the relationship between D/H and the average sightline gas density. 
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We use the simple model developed by Spitzer (1985) (see also Jenkins et al. 1986) to 
explain the correlation between the level of depletion and the mean line of sight density, to 
explore the relationship between D/H and (%)• In this model, the interstellar medium is 
composed of many warm diffuse clouds, higher density cold diffuse clouds which occur with 
a lower spatial frequency and cold dense clouds with even lower spatial frequency. At low 
densities, contributions from the diffuse warm gas dominates that of the other two types 
of clouds. Cold diffuse clouds start to prevail with increasing densities, and at even higher 
densities, contributions from cold dense clouds outweigh everything else. Therefore, any 
correlation of a physical property with (nn) can be attributed partly to differences of the 
properties of the types of clouds. 

Jenkins et al. (1986) used the model developed by Spitzer (1985) to derive Equation 
1, which predicts the abundance of element X i: A(X,j) = log N(Xi) — log N(R) (iV(H) = 
AT(HI) +2A^(H 2 )), as a function of the sightline properties: 

A(X t ) = A c (X t ) + log[l + -J^-Cf - 1)] (1) 

where S w (Xi) is the depletion of element X i in the warm gas and 8 c (Xi) in the cold clouds, 
and n w is the mean density of the warm gas along the line of sight 

Equation 1 can be simplified to 

A(X t ) = M*i) + M 1 + ^(10^^)-^)) - 1)] (2) 

where A w (Xi) and A c (Xi) are the abundances of Xi in the warm diffuse medium and cold 
clouds, respectively. This equation assumes that diffuse cold clouds and dense cold clouds 
have the same level of depletion and that for (%) / < 1, all the neutral hydrogen along 
the line of sight is in warm gas (note that for the gas densities considered here, cold dense 
clouds are not expected to contribute significantly to the sightline properties, Spitzer 1985). 
Hence, for (nn) < n w , A(Xi) has the constant value A w (Xi), characteristic of warm gas. 

Figure 19 (top panel) presents log(D/H) as a function of the average sightline density, 
(nn). For comparison, we plot also log(0/H) (bottom panel). All the data have been 
presented previously in the top and middle panels of Figure 17 and are also listed in Table 12. 
Also plotted is (D/H) pr i m = (2.62 ± o2o)xlO -5 ', determined from analysis of the Wilkinson 
Microwave Anisotropy Probe (WMAP) satellite data (top panel, Spergel et al. 2003) and 
O/H = (3.43 ± 0.15)xl0~ 4 from Meyer et al. (1998) (bottom panel). Like Jenkins et al. 
(1986) we adopt n w = 0.10 cm -3 below, to determine A W (D) and A C (D). As it can be seen 
in Figure 19, (jih) = 0.10 cm -3 seems to be a natural point, after which the scatter in D/H 
increases. However choosing a value in the range 0.08 < (n H ) < 0.12 does not affect A W (D) 
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determined below. A non-linear least squares fit (x 2 minimization) is used to fit A(D) from 
Equation 2 above to the data plotted in the top panel of Figure 19. 

We derive A W (D) = -4.80 ± 0.01 ((n H ) <0.10 cm' 3 ) corresponding to D/H = 1.58xl0 -5 
and xl — 1-2 (26 degrees of freedom). This D/H ratio is consistent with the one derived 
for the LB by Wood et al. (2004), D/H = (1.56 ± 0.04) xlO" 5 . We note however, that not 
all data points with (%) <0.10 cm~ 3 correspond to targets located inside the LB, although 
all sightlines inside the LB have (%) <0.10 cm" 3 . There are five sightlines of this type, 
outside the LB; (3 Cen, A Sco, 7 2 Vel, ( Pup, and TD1 32709. Three of these have D/H values 
that are consistent with the low density average but A Sco (D/HxlO 5 = 0.76 ± 0.25) and 
7 2 Vel (D/Hx 10 5 = 2.18 ± q 22 ,) are exceptional. Thus, for low densities, as well as for higher 
densities, which we discuss below, some sightlines apparently defy the trends by statistically 
significant margins. 

For (nn) >0.10 cm -3 , A C {D) is mostly affected by the five points with a high D/H 
ratio, i.e., log(D/H) > -4.8. We derive A C (D) = -5.03 (D/H = 0.93xl0~ 5 ) with xl = 4.7 
(16 degrees of freedom), using all ratios corresponding to sightlines with (n H ) >0.10 cm~ 3 
(solid line in Figure 19). If we assume that these five ratios are exceptional, corresponding 
to sightlines with different conditions, and remove them from our fit, we would derive A C (D) 
= —5.20 (D/H = 0.63xl0~ 5 ) with xl — 2.4 (11 degrees of freedom). The corresponding fit 
is represented by a dashed-dot line in the top panel of Figure 19. This fit is now heavily 
dominated by a single ratio corresponding to the BD+28 4211 sightline (with a small uncer- 
tainty compared to the other sightlines). If in addition we exclude this ratio so that the fit 
is not dominated by one point we would derive, A C (D) = —5.38 (D/H = 0.42 xlO -5 ) with 
xl = 1.5 (10 degrees of freedom), represented by a dashed line in the top panel of Figure 19. 

For O/H in the bottom panel of Figure 19 we derive the weighted mean of O/H = (3.78 
± 0.18) xlO -4 (1 a in the mean) for (%) < 0.1 cm~ 3 , with xl = 0-8 for 11 degrees of free- 
dom. The square-root of the weighted average variance (standard deviation) is 0.56 xl0~ 4 . 
Considering only the sightlines with (tir) > 0.1 cm -3 we derive the weighted mean O/H = 
(3.41 ± 0.15)xl0~ 4 with a standard deviation of 1.40 xl0~ 4 and xl = 5-8 for 14 degrees of 
freedom. The O/H data does not show the same trend with (ra H ) as the one displayed by 
D/H. Recent studies have shown that oxygen is not expected to be depleted for (%) < 1.5 
cm -3 (Cartledge et al. 2004). However, the bottom panel of Figure 19 does seem to indicate 
that for some of the densest sightlines there is an anti-correlation between D/H and O/H, 
which could be an indication of astration. 

Another interesting comparison is the depletion of Fe compared to that implied by 
Figure 19 for D. Jenkins et al. (1986) found that the abundance of Fe decreases by a factor 
of three from the warm gas value, by (nn) = 1-0. Figure 19 indicates a factor of two or more 



-31 - 



decrease over the same density range for D, supporting the correlation between the values 
of D/H and the Fe depletion reported by Linsky et al. (2005, in prep). 

Analyzing D/H as a function of the average sightline density removes the need for the 
three regimes displayed in Figure 17. In this picture, D/H is constant and has the same 
value as in the LB for (%) <0.10 cm -3 , i. e. for sightlines dominated by warm diffuse 
gas, implying that D /H is constant in a broader sense than has been considered before. For 
(tih) >0.10 cm' 3 , D/H decreases with increasing density, as more of the gas is associated 
with cold diffuse clouds. Considering that for other elements such as Si and Fe this has been 
considered as evidence of depletion of these elements into dust grains (see for e.g. Jenkins 
et al. 1986), then our study strongly suggests that D might also be depleted into dust grains 
as suggested by Draine (2004) (but see Linsky et al. 2005, in prep, for discussions on other 
possible explanations, such as variable astration and localized infall of metal-poor gas). 
However, as noted above, in both regimes there are a few sightlines which are exceptional, 
with high D/H ratios. It is possible that along these sightlines the conditions are such that 
D has either been released from the grains, does not deplete onto grains, or is only mildly 
depleted. Also, while there are no accurate quantitative estimates for significant production 
of D by local sources, the possibility of such sources cannot be entirely dismissed (see for 
example Mullan & Linsky 1999; Prodanovic & Fields 2003). The true value of D/H in 
cold diffuse clouds, A C (D), is difficult to derive at this point given the small number of 
sightlines with high gas density. Nevertheless, it seems reasonable to estimate that in cold 
diffuse clouds D/H < 0.93xl0~ 5 , which is the result of our fit when all the sightlines with 
(n H ) >0.10 cm" 3 are included in the analysis. 



8. SUMMARY 

We have used data obtained with FUSE together with archival data from IUE and 
HST to derive column densities and ratios of column densities along the lines of sight to 
WD 1034+001, BD+39 3226, and TD1 32709. The D/H derived here for two of these sight- 
lines are not consistent at the 1 a level with the Local Bubble value; none of the D/O ratios 
are consistent, at the 1 a level with the Local Bubble value, and present some of the highest 
values in the literature. Considered along with the ratios published for other sightlines, our 
results reinforce the scatter in the D/H measurements and indicate that there is also scatter 
in the D/O ratio, implying that the high D/H ratios derived here and along other sight- 
lines are unlikely to be due to problems with the A(HI) determinations. We estimate that 
the additional, presently unknown, systematic errors on the determination of A(HI) would 
have to be of the order of 40% of each D/H ratio, in order to bring all the D/H ratios into 
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agreement. Thus, our results support the idea that some sightlines in the Milky Way ISM 
have high D/H ratios. For sightlines with logA(H) > 20.7 the O/H ratios are not consistent 
with a single value, and are inconsistent with the local ISM value derived by Meyer et al. 
(1998). The hint of anti-correlation between D/H and O/H for these sightlines could be an 
indication of astration. 

The trend between decreasing D/H with decreasing Fe/H observed by Linsky et al. 
(2005, in prep) is also observed when D/O versus Fe/O are considered. These ratios are not 
subject to possible systematic uncertainties in the determination of N(H). 

In order to try to understand the behavior of the deuterium abundance we performed a 
study of D/H versus the sightline average gas density. Our results indicate that as long as the 
gas probed by the sightlines is in warm diffuse clouds with (%) <0.10 cm~ 3 D/H seems to 
be constant and have the Local Bubble value. In addition, our study shows that D/H seems 
to decrease with increasing sightline gas density, similarly to what has been observed for 
other elements such as Fe and Si, which strongly supports the idea that D might be depleted 
into dust grains. Finally, a few sightlines do not follow the trend, but show exceptionally 
high D/H ratios. 

This work is based on data obtained for the Guaranteed Time Team by the NASA- 
CNES-CSA FUSE mission operated by The Johns Hopkins University. Financial support 
to U. S. participants has been provided in part by NASA contract NAS5-32985 to Johns 
Hopkins University. Based on observations made with the NASA/ESA Hubble Space Tele- 
scope and the International Ultraviolet Explorer, obtained from the Data Archive at the 
Space Telescope Science Institute, which is operated by the Association of Universities for 
Research in Astronomy, Inc. under NASA contract NAS5-26555. Support for MAST for 
non-HST data is provided by the NASA Office of Space Science via grant NAG5-7584 and 
by other grants and contracts. The profile fitting procedure, Owens.f, used in this work was 
developed by M. Lemoine and the French FUSE Team. We thank Todd Tripp and Scott 
Friedman for their help in determining iV(H I) from IUE data and W. V. Dixon for useful 
discussions regarding the scattered light background in the FUSE data. We thank also Guil- 
laume Hebrard, Jeff Linsky, and Paule Sonnentrucker for comments that helped to improve 
the paper. 
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Table 1. Stellar Properties 



Star 


T e ff 


log g 


d 




I 


b 


Sp. Type 




— ^ISM 




(K) 


(cm -2 ) 


(pc) 




(°) 


(°) 




(km s x ) 


WD 1034+001 a 


100,000 ± iS:SSS 


7.5 ± 0.3 


155 ± 


58 
43 


247.55 


+47.75 


DO 




+50 


BD+39 3226 b 


45,000 


5.5 


290 ± 


140 
70 


65.00 


+28.77 


sdO 




-255 


TD132709 C 


46,500 ± 1,000 


5.55 ± 0.15 


520 ± 


90 


232.98 


+28.12 


sdO 




-13 d 



a From Werner et al. (1995). 
b From Bluhm et al. (1999). 
c From Dreizler (1993). 
d This work. 



Table 2. Log of FUSE observations 



Star 


Program ID 


Aperture 


Mode 


Time (Ks) 


Date 


CalFUSE 


WD 1034+001 


P1042003 


MDRS 


TTAG 


6.7 


2004 April 04 


2.4.2 


BD+39 3226 


P3020501 


MDRS 


HIST 


5.2 


2004 May 21 


2.4.2 


TD 132709 


P2051301 


MDRS 


HIST 


1.8 


2004 March 12 


2.4.1 




P2051302 


MDRS 


HIST 


3.2 


2004 May 4 


2.4.1 




P2051303 


MDRS 


HIST 


2.4 


2004 May 5 


2.4.1 
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Table 3. Atomic data and analysis methods for the lines used in the analyses" 



Species Wavelength (A) Log fX WD 1034+001 BD+39 3226 TD1 32709 



HI 


1215.670 


2.70 


P 




P 


DI 


916.1785 
917.8797 


-0.28 
-0.07 


P 




P 

A,P 




919.1013 


0.04 


A, C, P 


A, C, P 


A, CP 




920.713 


0.17 


A, C, P 


C,P 


A, CP 




922.899 


0.31 


A, C, P 




C 




925.974 


0.47 


C,P 


CP 


C, P 




930.495 


0.65 






C, P 




937.548 


0.86 


C 


c 


C, P 




949.485 


1.12 




c 


c 




972.272 


1.45 


c 


c 


c 


CII 


1036.3367 


2.09 


B/S 


A 


A 


CII* 


1037.018 


2.11 


A 


A 


A 


cm 


977.0200 


2.87 


A 


A 


B/S 


NI 


951.2948 


-1.66 






P 




952.303 


0.25 


C 


C 


C 




952.415 


0.21 




c 


c 




952.523 


-0.24 


A, C,P 


A, CP 


A, CP 




953.415 


1.10 


C 


C 


C 




953.655 


1.38 


c 


c 


c 




953.9699 


1.52 


c 


c 






954.1042 


0.81 


c 


c 


c 




955.8819 


-1.54 






p 




959.4937 


-1.30 






p 




963.990 


1.54 


c 


c 


c 




964.626 


0.96 


c 


c 


c 




1134.1653 


1.69 


c 


c 


c 




1134.4149 


1.53 


c 


c 


c 




1134.9803 


1.24 






c 


Nil 


1083.994 


2.10 


A 


A 


A 


01 


916.8150 


-0.36 


c 
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Table 3 — Continued 



Species Wavelength (A) Log fX WD 1034+001 BD+39 3226 TD1 32709 





919.658 


-0.06 


C 


c 


c 




919.917 


-0.79 


A,P 


A, CP 


A, C 




921.875 


0.04 


C 








922.200 


-0.65 


P 




p 




924.950 


0.15 


c 


c 






925.446 


-0.49 


A, C 


A, C 


C 




929.5168 


0.32 






C 




930.257 


-0.30 


C 


C 


C 




936.6295 


0.53 




c 


C 




948.6855 


0.77 






C 




950.885 


0.18 


C 


c 






971.738 


1.13 




c 


C 




974.070 


-1.82 


P, c 




P, c 




976.4481 


0.51 






c 




1039.2301 


0.98 






c 


OVI 


1031.9261 


2.13 


p 


p 


B/S 




1037.6167 


1.83 


p 




B/S 


Sill 


1020.699 


1.22 


A 


A 


A 


PII 


1152.8180 


2.45 


A 


A 


A 


S III 


1012.4950 


1.65 


A 






SIV 


1062.6780 


1.64 


A 






Arl 


1048.220 


2.44 


A 


A 


A 


Fell 


1055.2617 


0.90 




A, CP 






1063.176 


1.76 


C 


C 


C, P 




1081.9352 


-0.90 


C,P 


P 


C 




1096.8770 


1.55 


C,P 


C, P 


P 




1112.0480 


0.84 


CP 


p 


c 




1121.975 


1.36 


CP 


CP 






1125.448 


1.26 




CP 






1133.6654 


0.80 


A, CP 


CP 


A, C 




1142.3656 


0.68 


p 




C 
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Table 3 — Continued 



Species Wavelength (A) 


Log f\ 


WD 1034+001 


BD+39 3226 


TD 132709 


1143.2260 


1.31 


C,P 


CP 


A, CP 


1144.938 


2.08 


C,P 


C 





a A, C, and P, denote lines that are analyzed with apparent optical depth, curve 
of growth, and profile fitting techniques, respectively. B/S means that no column 
densities were derived from these lines due to blending with stellar features. 
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Table 4. Equivalent widths of the atomic lines used with the COG method (mA) 

Species Wavelength (A) WD 1034+001 BD+39 3226 TD1 32709 
DI 



NI 



01 



919.1013 


16 


96 


± 


2 


46 


12 


85 


± 


2 


72 


6.46 ± 4.23 


920.713 


18 


68 


± 


2 


44 


14 


68 


± 


3 


77 


18.59 


± 


2.23 


922.899 


27 


49 


± 


3 


47 












26.93 


± 


5.13 


925.974 


39 


78 


± 


3 


00 


26 


78 


± 


2 


60 


26.29 


± 


4.33 


930.495 






















41.78 


± 


3.10 


937.548 


46 


06 


± 


4 


42 


34 


77 


± 


3 


06 


48.31 


± 


6.60 


949.485 












53 


68 


± 


2 


88 


53.61 


± 


5.80 


972.272 


67 


64 


± 


8 


91 


71 


00 


± 


4 


91 


74.34 


± 


6.52 


952.303 


62 


29 


± 


3 


43 


37 


71 


± 


1 


68 


61.98 


± 


3.17 


952.415 












44 


11 


± 


1 


59 


49.23 


± 


2.52 


952.523 


19 


40 


± 


2 


71 


21 


96 


± 


1 


52 


29.63 


± 


2.82 


953.415 


67 


19 


± 


3 


58 


54 


64 


± 


2 


05 


70.61 


± 


3.48 


953.655 


75 


78 


± 


3 


84 


63 


32 


± 


2 


51 


81.25 


± 


3.53 


953.9699 


78 


27 


± 


3 


54 


65 


91 


± 


2 


14 








954.1042 


68 


76 


± 


3 


49 


51 


17 


± 


1 


86 


57.25 


± 


3.59 


963.990 


67 


68 


± 


4 


73 


77 


37 


± 


2 


54 


70.22 


± 


3.19 


964.626 


66 


84 


± 


4 


10 


64 


60 


± 


2 


50 


59.59 


± 


3.37 


1134.1653 


84 


50 


± 


2 


76 


73 


40 


± 


1 


31 


101.82 ± 2.35 


1134.4149 


86 


95 


± 


2 


48 


81 


27 


± 


1 


41 


107.53 ± 2.57 


1134.9803 






















97.89 


± 


2.01 


916.8150 


48 


33 


± 


3 


92 


















919.658 


56 


28 


± 


3 


32 


44 


98 


± 


2 


57 


44.90 


± 


3.31 


919.917 












18 


63 


± 


1 


47 


27.62 


± 


2.91 


921.875 


59 


46 


± 


3 


33 


















924.950 


66 


17 


± 


3 


59 


57 


68 


± 


2 


53 








925.446 


44 


46 


± 


3 


37 


39 


13 


± 


2 


14 


41.62 


± 


2.88 


929.5168 






















72.33 


± 


3.22 


930.257 


54 


78 


± 


3 


54 


41 


69 


± 


2 


11 


50.04 


± 


3.72 


936.6295 












70 


93 


± 


2 


69 


85.04 


± 


3.10 


948.6855 






















88.06 


± 


3.87 


950.885 


64 


30 


± 


3 


76 


54 


69 


± 


1 


75 
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Table 4 — Continued 



Species 


Wavelength (A) 


WD 1034+001 


BD+39 3226 


TD1 32709 




y /I. too 




y(J. (6 ± o.24 


OA OA _i_ k no 

84. o4 ± O.U2 




y ^4.u / u 


4.86 ± 2.U6 




o.yi ± i.yi 




n *7P. A A O 1 

y / 0.4481 






*7 A OO _l_ o o^; 




1 non oom 






1 n i n a _i_ i oo 
!Ul.y4 ± 1.88 


TTV^ TT 
re 11 


incc: 0£ 1 V 

lUoo.zol 1 




o.zy ± u.y4 






1063.176 


54.53 ± 2.81 


41.14 ± 1.28 


26.09 ± 1.23 




1081.9352 


21.03 ± 3.90 




8.00 ± 1.64 




1096.8770 


21.03 ± 1.75 


34.52 ± 1.43 






1112.0480 


11.19 ± 1.52 




5.87 ± 1.05 




1121.975 


22.58 ± 2.07 


14.14 ± 0.91 






1125.448 




22.41 ± 1.09 






1133.6654 


5.70 ± 1.90 


6.48 ± 0.64 


5.28 ± 1.01 




1142.3656 






5.84 ± 0.94 




1143.2260 


16.01 ± 1.92 


22.98 ± 1.06 


16.43 ± 1.47 




1144.938 


55.82 ± 2.48 


57.39 ± 1.34 





Table 5. H 2 Column Densities (Log) 



H 2 WD 1034+001 BD+39 3226 TD1 32709 



J = 





14.55 ± 0.10 


15.17 ± 0.15 


13.50 ± 0.15 


J = 


1 


15.63 ± 0.15 


15.41 ± °f 


14.35 ± 0.15 


J = 


2 


14.60 ± 0.10 


14.47 ± 0.10 


13.44 ± 0.15 


J = 


3 


14.30 ± 0.10 


13.95 ± 0.10 


13.33 ± 0.15 


J = 


4 




13.38 ± 0.10 




J = 


5 




13.49 ± 0.07 




Total 


15.72 ± gil 


15.65 ± 


14.48 ± ° il 
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Table 6. Column Densities (Log) Along the WD 1034+001 Sightline 



Species AOD COG PF 



DI 


15.27 ± ° i\ (A919) 


15.40 ± gig 


15.40 ± 


0.05 
0.03 


NI 


15.89 ± g;§g (A952.5) 


16.07 ± g;g 


15.84 ± 


0.12 
0.13 


01 


16.57 ± 0.07 (A919.9) 


16.62 ± 0.05 


16.57 ± 


0.07 
0.08 


Fell 


14.19 ± ° il (A1133) 


14 14 + °- n 


14.05 ± 


0.05 
0.03 



Note. - The transitions associated with the AOD results 
quoted in the table are given in (). Uncertainties are la. 



Table 7. Column Densities (Log) Along the BD+39 3226 Sightline 



Species AOD COG PF 



DI 


15.27 ± o 1 ;^ (A919) 


15.15 ± 


0.08 
0.07 


15.15 ± 8:^ 


NI 


15.93 ± 8:1? (A952.5) 


15.88 ± 


0.11 
0.09 


15.80 ± 8:o1 


OI 


16.33 ± 0.08 (A919.9) 


16.31 ± 


0.07 
0.06 


16.51 ± 0.04 


Fell 


14.05 ± 0.05 (A1055) 


14.13 ± 


0.06 
0.05 


14.21 ± 0.02 



Note. - The transitions associated with the AOD results 
quoted in the table are given in (). Uncertainties are la. 
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Table 8. Column Densities (Log) Along the TD1 32709 Sightline 



Species AOD COG PF 



DI 15.33 ± g;gg (A919) 15.29 ± ° ° l 15.26 ± 0.03 

NI 15.93 ± 0.04 (A952.5) 16.02 ± ° f 2 15.97 ± [];?? 

OI 16.36 ± g;g| (A919.9) 16.45 =b °;g§ 16.48 ± ggg 

Fell 13.98 ± (A1133) 14.03 ± 0.05 13.87 ± ° f 2 



Note. - The transitions associated with the AOD results 
quoted in the table are given in (). Uncertainties are la. 
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Table 9. Atomic Column Densities (Log) 



Species 


WD 1034+001 


BD+39 3226 


TD1 32709 


HI 


20.07 ± 0.07 


20.08 ± 0.09 a 


20.03 ± 0.10 


DI 


15.40 ± 0.07 


15.15 ± 0.05 


15.30 ± 0.05 


CII 


blend 


> 14.24 


> 14.65 


CII* 


> 13.88 


> 13.63 


> 13.94 


cm 


> 13.66 


> 13.40 


blend 


NI 


15.96 ± 0.12 


15.85 ± 0.10 


15.98 ± 0.10 


TVT TT 

IN 11 


S 1/( u 

^ 14.10 


14. AO 


14. 60 


01 


16.60 ± 0.10 


16.40 ± 0.10 


16.42 ± 0.10 


OVI 


13.08 ± 0.08 c 


13.14 ± 0.05 




Sill 


> 14.68 


> 14.58 


> 14.66 


PII 


> 13.16 


> 13.01 


> 13.19 


S III 


> 14.17 c 






SIV 


> 14.08 c 






Arl 


> 13.65 


> 13.45 


> 13.64 


Fell 


14.10 ± 0.10 


14.15 ± 0.07 


13.95 ± 0.10 



a From Bluhm et al. (1999). 

b Possibly contaminated by stellar Nil absorption. 
c Circumstellar contribution possible, see §4.2.1. 



Table 10. Ratios of Column Densities 



Ratio WD 1034+001 BD+39 3226 TD1 32709 



(DI/HI)xl0 5 


2.14 


± 


0.53 
0.45 


1.17 


± 


0.31 
0.25 


1.86 


± 


0.53 
0.43 


(NI/HI)xl0 5 


7.76 


± 


2.82 
2.20 


5.89 


± 


2.04 
1.63 


8.91 


± 


3.26 
2.59 


(OI/H I)xl0 4 


3.39 


± 


1.06 
0.86 


2.09 


± 


0.72 
0.58 


2.45 


± 


0.90 
0.71 


(DI/OI)xl0 2 


6.31 


± 


1.79 
1.38 


5.62 


± 


1.61 
1.31 


7.59 


± 


2.17 
1.76 


(DI/NljxlO 1 


2.75 


± 


1.00 
0.78 


2.00 


± 


0.57 
0.46 


2.09 


± 


0.60 
0.49 


OI/NI 


4.37 


± 


1.79 
1.38 


3.55 


± 


1.30 
1.03 


2.75 


± 


1.01 
0.80 



-47- 



Table 11. BD+39 3226: N FUSE vs. N rfeus+we 



Species log N FUSE a log N ORFEUS+IUE h 



DI 






15.15 ± 0.05 


15.16 ± 0.13 


C 






> 14.38 


16.40 ± 0.75 


NI 






15.85 ± 0.10 


14.75 ± 0.25 


01 






16.40 ± 0.10 


16.40 ± g : ™ 


Si 






> 14.58 


14.80 ± 0.20 


Fell 




14.15 ± 0.07 


14.10 ± 0.15 


H 2 


{J 


= 0) 


15.17 ± 0.15 


15.10 ± 0.20 




{J 


= 1) 


15.41 ± Ho 


15.50 ± 8:18 




(J 


= 2) 


14.47 ± 0.10 


i4.io ± 




(J 


= 3) 


13.95 ± 0.10 


13.95 ± ; 10 




(J 


= 4) 


13.38 ± 0.10 


13.80 ± 8:18 




(J 


= 5) 


13.49 ± 0.07 


< 13.40 



a This work. 

b Bluhm et al. (1999). No /-value corrections 
were applied to N. 



Table 12. Compilation of Sightline Ratios 



Star d(pc) logiV(H) (D/H)xl0 5 (D/O)xl0 2 (O/H)xl0 4 (n H > (cm" 3 ) Ref. 



Sirius B 




2.64 


17.81 


± 


0.11 


1.17 ± 


0.37 


3.9 ± 


0.8 


3.02 ± 


0.95 
0.75 


0.079 


2, 2,1 


e Eri 




3.22 


17.88 


± 


0.04 


1.40 ± 


0.20 












0.076 


3, 3,- 


Procyon 




3.50 


18.06 


± 


0.05 


1.60 ± 


0.20 












0.106 


4 


e Ind 




3.63 


18.00 


± 


0.05 


1.60 ± 


0.20 












0.089 


5 


36 Oph 




5.99 


17.85 


± 


0.08 


1.50 ± 


0.25 












0.038 


6 


P Gem 




10.34 


18.26 


± 


0.04 


1.47 ± 


0.20 












0.057 


3, 3,- 


Capella 




12.9 


18.24 


± 


0.04 


1.60 ± 


0.14 
0.19 


2.6 ± 


1 


2 


6.31 ± 


3.69 
1.95 


0.044 


4, 4,7 


P Cas 




16.7 


18.13 


± 


0.03 


1.70 ± 


0.15 












0.026 


3 


a Tri 




19.7 


18.33 


± 


0.04 


1.32 ± 


0.30 












0.035 


3 


A And 




25.8 


18.45 


± 


0.08 


1.70 ± 


0.25 












0.035 


7 


P Cet 




29.4 


18.36 


± 


0.05 


2.20 ± 


0.55 












0.025 


8 


HR1099 




29 


18.13 


± 


0.02 


1.46 ± 


0.09 












0.015 


8 


cr Gem 




37 


18.20 


± 


0.04 


1.36 ± 


0.20 












0.014 


3 


WD 1634- 


-573 


37 


18.85 


± 


0.06 


1.60 ± 


0.25 


3.5 ± 





3 


4.57 ± 0.75 


0.062 


9 


WD2211- 


-495 


53 


18.76 


± 


0.15 


1.51 ± 


0.60 


4.0 ± 





6 


3.80 ± 


1.61 
1.16 


0.035 


10 


HZ43 




68 


17.93 


± 


0.03 


1.66 ± 


0.14 


4.6 ± 





5 


3.63 ± 


0.44 
0.70 


0.004 


11 


G191-B2B 


69 


18.18 


± 


0.09 


1.66 ± 


0.45 


3.5 ± 





4 


4.79 ± 


1.19 
0.99 


0.007 


12 


Feige 24 




74 


18.47 


± 


0.03 


1.30 ± 


0.50 


3.9 ± 


1 


7 


3.33 ± 


0.59 
0.56 


0.013 


13, 13, 14 


WD 0621- 


-376 


78 


18.70 


± 


0.15 


1.41 ± 


0.56 


3.9 ± 





6 


3.63 ± 


1.56 
1.10 


0.021 


15 


GD246 




79 


19.11 


± 


0.03 


1.51 ± 


0.20 
0.17 


4.2 ± 





6 


3.63 ± 


0.44 
0.34 


0.053 


16 


a Vir 




80 


19.00 


± 


0.10 


1.58 ± 


1.01 
0.46 


4.2 ± 


2 





3.80 ± 


1.39 
1.11 


0.041 


17, 17, 18 


31 Com 




94 


17.88 


± 


0.03 


2.00 ± 


0.20 












0.003 


3 



Table 12— Continued 



Star 


d (pc) 


log AT(H) 


(D/H)xl0 5 (D/O)xl0 2 (O/H)xl0 4 ( 


n H ) (cm 3 


) Ref. 


a Cru 


no 
yo 


1 n p.c\ _i_ n 1 n 

ly.ou ± U.1U 


2.24 ± ° e 4 

u. oz 






n 1 qo 
U.loZ 


1 7 1*7 
1 / , 1 f , - 


nn i oc /loi 1 
JdU+zo 4Z11 


i n a 
1U4 


m or: _i_ n no 

ly.oo ± u.uz 


1.26 ± 0.08 


4.7 ± 0.4 


2.69 ZTZ n oo 


n oo 1 
U.ZZ1 


19, 1, 1 


y oar 


loo 


on os i n 1 n 

ZU.Zo ± U.1U 


0.50 ± 0.16 






n /I K7 
U.40 l 


on 

ZU 


p ima 


lOo 


i o on i n 1 p. 
lo.ZU ± U.lo 


!-20 ± li° 

U.OU 






n nnQ 
U.UUo 


01 
Zl 


"wr\ i nQ/i i nm 
WJJ !Uo4+UUl 


100 


on n7 -i- n n7 
ZU.U / ± U.U / 


2.14 ± ° f 5 

U.40 


6.31 ± \ H 


3.39 ZTZ ri oc 

U.oO 


n o/i f; 
U.Z40 


oo 
zz 


p Uen 


lol 


in e.A -i- n n k 
iy.04 ± U.UO 


1.26 ± PI 

u.4o 






n n7n 
U.U /U 


1 7 
1 / 


reige liu 


i 70 

i ( y 


on 1 /i -i- n nn 

zu.i4 ± u.uy 


2.14 ± 0.41 


2.6 ± 0.5 


8.32 ± ?'qj 


n o/i k 
U.Z40 


OQ OQ OK 

Zo, Zo, oO 


7 Cas 


loo 


on n/i -i- n n/i 
ZU.U4 ± U.U4 


1.12 ± 0.25 


2.5 ± 0.4 


5.25 i q ^ 


n i qo 

u. toy 


O/I O/I OK 

Z4, Z4, ZO 


A OCO 


01 
Z10 


i n o'? -i- n oq 
iy. zo it u.Uo 


0.76 ± 0.25 


1.8 ± 81 


4.24 ± 81? 


n nofi 

U.UZO 


Ofi 
ZD 


7 vei 


ZOo 


i o 71 _i_ n oq 
iy. I i ± U.Uo 


2.18 ± ° H 






n nR/i 
U.U04 


07 07 

ZY, ZY, - 


o uri 


081 
Z01 


on 1 q -i- n n°. 
zu.iy it u.uo 


0.74 ± 

u.uy 


2.5 ± 0.4 


2.80 ± 0.40 


n 1 7o 
u. i / y 


08 os ok. 
zo, zo, zo 




oon 
zyu 


on ns -U n no 
zu.Uo it u.uy 


1.17 ± °i 

U.zo 


5.62 ± } " 


2.09 ± ° H 

U.OO 


n i Q/i 

U. 104 


oo oo oo 
zy, zz, zz 


tdo nnoo i 1 on 


on7 
zy / 


on aq _i_ n n /i 
ZU.4o it U.U4 


1.91 ± g;» 


2-4 ± H 


7-76 ± lf 3 


n QQn 
U.ooU 


on 

oU 


t Uri 


4U / 


on 1 £ _i_ n 1 n 
ZU.lO it U.1U 


1 40 ± 9-S 
- L - 4:U ^ 1.00 


3.5 ± 0.8 


3 98 ± H? 

kj • %J \J _i — i ii 


U.llO 


Q 1 Q 1 OK 

ol, ol, Zo 




41 9 


on An _i_ n no 
z,u.^±u in u.uo 


0.65 ± 0.30 


1.9 ± 0.3 


3.80 ± 8:?8 


u. lyo 


O J. , O -L , 


C Pup 


429 


19.96 ± 0.03 


1.42 ± °;* 






0.069 


27 


TD 132709 


520 


20.03 ± 0.10 


i.86 ± 8:1 


7.59 ± ?;« 


2.45 ± 8'?? 


0.067 


22 


LSE44 


554 


20.52 ± ° \l 


2 24 + °- 70 


i.99 ± 8:1 


ii-3 ± a 


0.194 


32 


LS 1274 


580 


20.98 ± 0.04 


0.76 ± 0.18 


1.8 ± 0.5 


4.68 ± HI 


0.534 


33, 33, 1 


JL9 


590 


20.78 ± 0.05 


1.00 ± 0.19 


1.9 ± 0.8 


5.25 i 


0.331 


33 


HD 195965 


794 


20.95 ± 0.03 


0.85 ± g;g 


1.3 ± 0.3 


6 61+ °- 52 
u.ui it 56 


0.364 


34 


HD 191877 


2200 


21.05 ± 0.05 


o.78 ± 8:3 


2.5 ± 1.0 


q nn _l 0.99 
o.uy it 049 


0.165 


34 


HD 90087 


2740 


21.17 ± 0.05 


0.98 ± 0.19 


1.7 ± 0.4 


5.8 ± 1.0 


0.175 


35 



Table 12— Continued 



Star d(pc) logiV(H) (D/H)xl0 5 (D/O)xl0 2 (O/H)xl0 4 (n H ) (cm~ 3 ) Ref. 



WD 2004-605 


52 


19.11 


4.8 ± 1.3 


1 


WD 1631+781 


67 


19.36 


3.8 ± 1.2 


1 


WD 2331-475 


82 


18.94 


5.1 ± 1.1 


16 


HZ 21 


115 


19.20 


4.6 ± 1.0 


16 


Lan23 


122 


20.18 


3.2 ± 1.6 


16 


CPD-31 1701 


131 


19.39 


3.7 ±0.6 


1 



Note. - The reference numbers refer to iV(H), N(D), and N(0), respectively. When 
only one reference is given it refers to all the N quoted for that particular sightline. The 
second part of the table contains sightlines for which only D/O is available. These ratios are 
displayed in Figure 17 for comparison purposes only. For these sightlines the quoted AT(H) 
are estimated by using N(0) and O/H = (3.43 ± 0.15)xl0~ 4 from Meyer et al. (1998). 

References. - (1) Hebrard & Moos (2003), (2) Hebrard et al. (1999), (3) Dring et al. 
(1997), (4) Linsky et al. (1995), (5) Wood et al. (1996), (6) Wood et al. (2000), (7) Wood 
et al. (2002b), (8) Piskunov et al. (1997), (9) Wood et al. (2002a), (10) Hebrard et al. (2002), 
(11) Kruk et al. (2002), (12) Lemoine et al. (2002), (13) Vennes et al. (2000), (14) Oliveira 
et al. (2005), (15) Lehner et al. (2002), (16) Oliveira et al. (2003), (17) York & Rogerson 
(1976), (18) York & Kinahan (1979), (19) Sonneborn et al. (2002), (20) Allen et al. (1992), 
(21) Gry et al. (1985), (22) This work, (23) Friedman et al. (2002), (24) Ferlet et al. (1980), 



(25) Meyer et al. (1998), (26) York (1983), (27) Sonneborn et al. (2000), (28) Jenkins et al. 
(1999), (29) Bluhm et al. (1999), (30) Williger et al. (2005), (31) Laurent et al. (1979), (32) 
Friedman et al. (2005), (33) Wood et al. (2004), (34) Hoopes et al. (2003), (35) Hebrard 
et al. (2005). 
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Fig. 1. — FUSE spectrum of WD 1034+001 with identification of interstellar and stellar 
absorptions. Labels in [] refer to stellar lines; dashed vertical lines mark the position of the 
strongest H 2 lines along this sightline. Stellar HI absorption is displaced from the ISM one 
by the velocity shift specified in Table 1. Note the isolated broad Hell stellar lines around 
992 and 958 A. Other stellar He II lines fall close to the positions of the Lya series and are 
not resolved. The FUSE channel used for each panel is indicated at the bottom right. The 
spectrum was binned by 4 for display purposes only. 
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Same as Fig. 1 but for the BD+39 3226 3 sightlme / 
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Fig. 3.— Same as Fig. 1 but for the TD1 32709 sightline. 
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Fig. 4. — H 2 excitation diagram for the three sightlines. Each J level is labeled at the top 
of the first panel. See §4.1 for discussion. 
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Fig. 5. — Curve of growth for 01 along the WD 1034+001 sightline. We derive log iV(OI) 
= 16.62 ± 0.05. Dashed lines indicate the 1 a uncertainty in N. 
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Fig. 6. — Fits to some of the lines used in the analysis of the WD 1034+001 sightline. The 
FUSE channels corresponding to the plotted data are indicated in the bottom left of each 
panel. For H 2 , the corresponding J level is indicated in parenthesis. [] are used to indicate 
lines of photospheric origin. 
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Fig. 7. — Fell A1144.9 interstellar absorption toward BD+39 3226 as seen with FUSE. Only 
one absorption component, at v® = —25 km s^ 1 , is present. We see no indication of a weaker 
absorption component at v & = —75 km s _1 , as reported by Bluhm et al. (1999). 
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Fig. 8.— Same as Fig. 6 but for the BD+39 3226 sightline. 
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Fig. 9. — Curve-of-growth for O I along the BD+39 3226 sightline with the la uncertainties 
in TV plotted as dashed lines. The fit yields log iV(OI) = 16.31 ± g;g|. 
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Fig. 10. — Comparison between the SiC IB (solid) and SiC 2A (dotted) data for the 
BD+39 3226 sightline in the region of the OI A974 transition. Data at the top were cal- 
ibrated with version 3 of CalFUSE (V3), data at the bottom were calibrated with version 
2.4 (V2.4) (data shifted vertically for clarity). The apparent discrepancy between the profiles 
around O I A974 in the V2.4 SiC IB and SiC 2A data disappears when V3 is used. See §4.2.2 
for discussion. 
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Fig. 11.— Curve of growth for 01 along the TD1 32709 sightline. We derive log iV(OI) = 
16.45 ± [];q3. Dashed lines indicate the 1 a uncertainty in N. 
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Fig. 12.— Same as Fig. 6 but for the TD1 32709 sightline. 
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Fig. 13. — Lya GHRS observations of the WD 1034+001 sightline. Three stellar models are 
overplotted: T e g = 100,000 K, logg = 7.5 (best fit model, solid line), T c g = 115,000 K, logg = 
7.2 (dotted line), and T c r = 90,000 K, logg = 7.8 (dash-dotted line). All the models have 
log (He/H) = 3.0. The two last models produce the highest and lowest iV(HI), respectively, 
and are used to determine the uncertainties in AT(HI) associated with the stellar models. 
Geocoronal emission is annotated with ©. We derive log iV(HI) = 20.07 ± 0.07 (lcr). 
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Fig. 14. — Fit to the Lya interstellar absorption along the WD 1034+001 sightline, using 
the best fit stellar model to normalize the data. Terrestrial airglow is labeled with ©. 
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Fig. 15. — Lya IUE observations of the TD1 32709 sightline. Three stellar models are 
overplotted: T e g- = 46,500 K, \ogg = 5.55 (best fit model, solid line), T e g = 45,000 K, logg = 
5.7 (dotted line), and T e g = 47500 K, logg = 5.4 (dash-dotted line). All the models have 
log (He/H) = 2.0. The two last models produce the lowest and highest iV(HI), respectively, 
and are used to determine the uncertainties in AT(HI) associated with the stellar models. 
Geocoronal emission is annotated with ©. We derive log iV(HT) = 20.03 ± 0.10 (lcr). 
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Fig. 16. — Fit to the Lya interstellar absorption along the TD1 32709 sightline, using the 
best fit stellar model to normalize the data. Top panel: The blue wing of the Lya line is 
used to define the zero flux level, yielding logiV(HI) = 20.03. Bottom panel: The red wing 
of the Lya line is used to define the zero flux level, yielding logiV(HI) = 20.12, within the 1 
a uncertainty quoted in Table 9 for iV(HT) along this sightline. 
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Fig. 17. — D/H, 0/H, and D/0 ratios as a function of the sightline total hydrogen column 
density, iV(H). Sightlines for which all three D/H, 0/H, and D/0 ratios are available are 
marked by asterisks (literature values) and filled circles (this work). Top panel: Sightlines 
for which no 1 measurement is available are marked with squares, the error bars are 
displayed by dotted lines. (D/H) lb = (1.56± 0.04) xl0~ 5 (derived by Wood et al. 2004, from 
a compilation of values in the literature) is indicated by solid and dashed (la) horizontal 
lines for log iV(H) < 19.2, which defines the contour of the Local Bubble (dashed vertical 
line). For 19.2 < log iV(H) < 20.7 there is a large scatter in the D/H ratio. For log iV(H) 
> 20.7, D/H seems to be constant, but lower than the LB value. Middle panel: The O/H 
= (3.43 ± 0.15) xl0~ 4 derived by Meyer et al. (1998) is marked by solid and dashed (la) 
horizontal lines. Bottom panel: D/O = (3.84 ± 0.16)xl0~ 2 derived by Hebrard & Moos 
(2003) for the LB is marked by solid and dashed (la) horizontal lines. Sightlines for which 
no N(H) is available are indicated by squares, the uncertainties are displayed by dotted 
lines. For these sightlines we use iV(OI) and the O/H ratio derived by Meyer et al. (1998) 
to estimate N(H). 
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Fig. 18. — D/0 as a function of Fe/O. The solid line represents the fit to all the data points 
presented. The dashed line represents the fit to the data when the four data points displayed 
by squares are removed from the fit. See §7.2 for discussion. 
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Fig. 19. — Top panel: Log(D/H) as a function of the average sightline density, (%) (cm -3 ). 
The WM^P-based result of (D/H) prim = (2.62 ± [J;^)xl0~ 5 (Spergel et al. 2003) is also 
plotted (solid and dotted lines at log(D/H) = —4.57). The solid line represents the fit of 
Equation 2 to all data points, yielding A W (D) = -4.80 (D/H = 1.58xl0 -5 ) and A C (D) = 
-5.03 (D/H = 0.93xl0 -5 ). The dash- dotted and dashed lines correspond to A C (D) = -5.20 
(D/H = 0.63 xl0~ 5 ) and A C (D) = -5.38 (D/H = 0.42 xlO" 5 ), respectively. See §7.3 for 
discussion. Bottom panel: Log(0/H) as a function of (%). The solid and dashed lines 
represent O/H = (3.43 ± 0.15) x!0~ 4 derived by Meyer et al. (1998). 



